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Abstract 


This report describes a viable, efficient, and easily mechanized carrier regenerat- 
ing receiver for use in suppressed carrier-tracking systems. The receiver referred 
to as a data-aided receiver (DAR) incorporates a data-aided loop (DAL) which 
provides the required carrier reference signal. The DAL employs the principle of 
decision feedback and as such is more efficient than other forms of suppressed 
carrier-tracking loops 

The analysis, design, and implementation of the DAR are covered in detail Per- 
formance comparisons and mechanization tradeoffs are made, wherever possible, 
with discrete carrier systems and other suppressed carrier systems presently in 
use. Experimental performance verification is given throughout m support of the 
theory presented. 
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The Theory, Design, and Operation of the Suppressed 
Carrier Data-Aided Tracking Receiver 


I. Introduction 


Time and progress have made tremendous inroads on 
the space-oriented communication system technology that 
existed in 1958. Data bit rates have increased some four 
orders of magnitude, and communication distances pres- 
ently reach to the outer planets. Ground antennas have 
grown from the relatively small helical types used for the 
early “Microlock” system to giant 64-m parabolic struc- 
tures with gains of 72 dB at X-band Receiver noise tem- 
peratures have dropped from over 1000 K to less than 25 
Data processors have evolved from hand reduction based 
on strip-chart records to large computer systems. 

This progress, of course, took time and research. Ideas 
had to be generated and developments made until the 
technology reached fruition. Often, the implementation 
and acceptance of a new advancement has been paced by 
the availability of a good, reliable, low-cost circuit tech- 
nology. The use of block and convolutional codes certainly 
fits this pattern as their properties and advantages were 
known well before their actual use beginning in the late 
’60s, when integrated circuits and minicomputers made 
decoding practical. 

The cost of further improvements in present communi- 
cations capability, in general, appears to be very high. 
A great amount is invested in equipment, so much so that 


to do other than make small modifications within the exist- 
ing systems is economically unfeasible. With this as a 
premise, emphasis is placed on ways by which the system 
performance can be improved without major modification 
or replacement. 

One such area for improvement involves the amount of 
transmitter power allocated to the carrier. From the incep- 
tion of the U.S. space program in the late 1950s, the use 
of carrier coherent receivers has formed the basis for 
tracking and communicating with the space vehicle At 
the heart of the receiver is the phase-locked loop (PLL), 
which even today continues to be constructed and used 
almost identically to that used for the first Explorer 
satellite 

The concept of using what is called “suppressed carrier” 
signaling to improve overall system performance, as com- 
pared to “discrete carrier” signaling, is very old The 
advantages of such an approach as it pertains to modem 
digital communications are presented in Section II of this 
report 

The acceptance of suppressed carrier signals for space 
vehicles has generally been thwarted by the lack of a com- 
plete cohesive theory which explains how well such a 
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signaling system, including its “carrier regenerating re- 
ceiver,” works relative to the well established discrete 
carrier systems. The use of carrier regenerating receivers 
has been further delayed due to the lack of a viable, effi- 
cient, and easily mechanized design. This report, which 
contains complete analysis, designs and implementations, 
and applications for a receiver known as the data-aided 
receiver (DAR), seeks to eliminate the above deficiencies 

Section III of the report deals with the general theory 
of operation for the suppressed carrier DAR, in both the 
acquisition and steady-state tracking modes Functions 
and forms which are synonymous with the PLL are fully 
discussed. 

Section IV concentrates on the detailed implementation 
of the DAR, with mechanization tradeoffs, analysis, and 
experimental performance verification. Section V, which 
is closely allied with Section IV, investigates the auto- 


matic gain control (AGC) and lock detection aspects of 
integrating the DAR into a larger receiver system. 

Section VI explores a number of practical applications 
and configurations of the DAR. Included are ways of 
modifying discrete carrier receivers to become data-aided 
receivers, applications to telemetry and command links, 
including subcarrier demodulation, and two-way tracking 
systems. 

The final section, VII, presents a qualitative summary 
of the advantages and disadvantages of DAR systems. 

Although no report can be absolutely complete, answer 
all questions, or cover all possible applications, the authors 
sincerely believe that this document represents a solid 
basis upon which one can decide whether or not the use 
of suppressed carrier signaling is both practical and eco- 
nomical for a given application 
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II. Advantages of Suppressed Carrier Signaling 


A. Basic Rationale— Power Division and 
Bandwidth 

For space applications, two major functions are gen- 
erally supplied by the radio link: (1) navigation via coher- 
ent tracking and (2) information in the form of telemetry 
or commands. In nearly all past programs, the tracking has 
been with respect to a discrete carrier component of the 
transmitter signal, and telemetry (or command) has been 
supplied by biphase-modulating data onto one or two sub- 
carriers, which in turn phase-modulate the transmitter. 
Expressed mathematically, for the single subcarrier case, 
the transmitter signal is 

s ( t ) = \T2 P t cos {u c t + 6 m d ( t ) 5 (f)} 

= y[2P T cos 0 m cos <o c t — \ 2 P T d ( t ) 5 (t) sin Q n sin « c t 
= V2F 0 cos a c t — V 2 P D d ( t ) 5 (t) sin o> c t (2-1) 

where 


Sections 5.4.5, 5 6, and A5.1). Let us define the efficiency 
of the discrete carrier system as 

t ?dc = -^ = sin 2 6 m (2-2) 

Ct 

As 7r/2, P c -> 0 and 77 DC -> 1. When 6 m = tt/2, P c = 0 
and there is no discrete carrier component, hence the term 
suppressed carrier. Now all of the available transmitter 
power is in the information-bearing signal, the best of 
conditions from the data signal point of view. If we define 
the efficiency of the suppressed carrier system rj sc also as 
the ratio of data signal power to total signal power, then 
clearly rise — 1. Finally, we may define the data signal 
improvement of the suppressed carrier format relative to 
that of the discrete carrier, with all other parameters held 
constant, as 

C D tm = CSC 2 Q m (2.3) 

Vac 


P T — total transmitter power 

P c = Pt cos 2 6 m — discrete carrier power 

P D = Pt sin 2 0 m = data signal power 

mc = carrier angular frequency 

d ( t ) = data stream with ±1 symbols 

5(f) = 5 in (aset) = unit power square- wave 
subcarrier of angular 
frequency « > s r 

Square-wave rather than sine-wave subcarriers are gen- 
erally used because of their higher efficiency with respect 
to sideband power utilization and ease of implementation 
on the spacecraft, thus, in the following discussion, square- 
wave subcarriers are assumed. From Eq. (2-1), it can be 
noted that the transmitter signal consists of two distinct 
parts: the discrete carrier (left-hand term), and the data 
signal (right-hand term), which is also referred to as the 
sideband term. From the information signaling point of 
view, the power devoted to the discrete carrier term is 
wasted. However, conventional discrete carrier-tracking 
receivers require a certain amount of carrier power in 
order to satisfactorily accomplish the navigation and 
carrier demodulation functions In fact, there exists a 
power division which optimizes the performance of both 
the tracking and data detection functions (see Ref. 2-1, 


Noting that 1 < G n < 00 , some data signal improve- 
ment is always achieved by using the suppressed carrier 
technique. 

With the suppressed carrier signal, the need to “carrier 
track” still exists. Since no discrete carrier signal is avail- 
able, however, a receiver which tracks the information- 
bearing signal is required In general terms this is called 
a carrier regenerating loop Several ways in which such a 
loop may be accomplished are introduced in Section C 
At this point, let us confine the discussion to carrier-track- 
ing efficiencies A figure of merit for the carrier-tracking 
loop is the signal-to-noise ratio (SNR) p in the carrier- 
tracking loop noise bandwidth W t For the discrete 
carrier-tracking loop, this signal-to-noise ratio is given by 

2 P c _ 2 P t cos 2 6 m 

Pnc ~ N 0 W,. - N n \V r , (2 ' 4) 

where N 0 is the one-sided noise spectral density in watts 
per hertz (W/Hz). For the carrier regeneration tracking 
loop, the loop SNR is given by 


2 P tLr 

psc = n 0 w l 


(2-5) 


where L, c < 1 is an efficiency loss associated with the 
carrier regeneration process Loss L R can range from 
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0.99 + for data-aided loops to less than 0.2 for squaring- 
type loops (see Section C). We may therefore define the 
tracking loop improvement of the suppressed carrier rela- 
tive to the discrete carrier case as 

G TR = — = L R sec 2 6 m (2-6) 

Pdc 

To give a feeling for the magnitude of the data signal 
and tracking loop improvements, Table 2-1 shows some 
typical values. ( L R = 0.91 is assumed.) 


Table 2-1. Typical values for data signal and tracking 
loop gains (Ls = 0.91) 


e m , deg 

Go, dB 

Gtr, dB 

30 

60 

08 

40 

38 

18 

50 

22 

34 

60 

1.2 

60 

70 

06 

94 


The foregoing discussion is, of course, somewhat simpli- 
fied, as other considerations and tradeoffs are necessary; 
however, it does give a gross feeling for what can be 
accomplished with suppressed carrier signals 

The suppressed carrier concept can easily be extended 
to the case of two data-modulating signals. The transmitter 
signal for the interplex* form of modulation (Ref. 2-1, Sec- 
tion 5.4.3) is given as 

s (£) = V 2fV sin [<o c t + 6 m idi (t) (t) 

+ e m2 d 1 (t)d 2 (t)S 1 (t)S 2 (t)] (2-7) 

The power in each usable component is given in Ref. 2-1 
(Table 5-6). When 6 ml — it/ 2, the suppressed carrier case 
is obtained, so that 

s ( t ) — V 2 P T cos 0 m2 di (t)5i (t) cos 

+ V 2 P T sin 8 m2 d 2 ( t ) 5 2 (t) sin m c t (2-8) 

It is seen that there are only two data-dependent compo- 
nents. A special, well known case, referred to as quadri- 
phase, is obtained when 6 m2 = tt/4, 5j ( t ) = 5 2 ( t ) = 1, and 
di (#) and d 2 (t) are derived from a single data stream. 


*This form is necessary to obtain quadrature suppressed earners 


Gains identical to those of the single-channel case are 
easily obtained as 

Gfli = = esc 2 d m i (2-9) 

“o 1 

G D2 = - = CSC 2 e ml (2-10) 

* 1)2 

G tr = = L/t sec 2 0 ml (2-11) 

* DC 

Equation (2-11) assumes that the carrier regeneration loop 
tracks the signal component dependent on d, (f). 

We now turn our attention to spectral occupancy. In 
the case of telemetry, the discrete carrier signal must be 
compatible with the tracking stations. For deep space 
applications interfacing with the Deep Space Network 
(DSN), there are restrictions on the choice of subcarrier 
frequency and the data symbol rate 


An unwritten but well observed restriction is that the 
minimum number of square-wave subcarrier cycles per 
symbol be 1.5. Actually, 3 or more subcarrier cycles per 
symbol is preferred. Another requirement is that the 
receiving system be able to recover up to and including 
the ninth harmonic of the subcarrier in order to maintain 
reasonable detection efficiency. A typical single-channel 
power spectrum, S («>), for the case of 1.5 subcarrier cycles 
per symbol is shown in Fig. 2-1. The radio-frequency (RF) 
bandwidth required to recover the sideband power inclu- 
sive of the ninth harmonic of the subcarrier is approxi- 
mately 20 o>sc — 30(2ir/R s ). 

For the suppressed carrier, single-channel signal, a sub- 
carrier is not required and, in fact, is a detriment to the 
carrier regenerating loop operation. Thus, the carrier may 
be directly modulated by the symbols; a typical spectrum 



Fig. 2-1. Typical discrete carrier single-channel spectrum 
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Fig. 2-2. Typical suppressed carrier single-channel spectrum 

is shown in Fig. 2-2. For this spectrum, the amount of RF 
bandwidth required to recover the same amount of side- 
band power as that from the discrete carrier spectrum of 
Fig. 2-1 is approximately 5(2: t£R s ). As a result, an RF 
bandwidth saving of a factor of 6 is obtained when the 
suppressed carrier signal is used. 

Similar bandwidth reduction arguments can be made 
for the two-channel case In general, the highest data- 
rate stream directly modulates the carrier. However, 
since there are several useful variations involved in the 
choice of modulation formats, the details will be omitted 
here. Section VII provides some insight into the possible 
variations. 

The general conclusion reached from the foregoing dis- 
cussion is that suppressed carrier modulation is highly 
efficient from the standpoint of both transmitter power 
utilization and RF bandwidth required. 

B. Two-Way SNR Advantages Relative to 
Discrete Carrier-Tracking Receivers 

At this point, we will briefly consider the merits of using 
suppressed carrier receivers (SCRs) in the two-way track- 
ing system. The tracking and data channel gains discussed 
in the preceding section apply to both up- and downlinks 
Our concern here is with the turnaround function. In the 
standard discrete carrier transponder, the uplink carrier is 
tracked by a phase-locked loop (PLL) (Ref. 2-2, Sections 
2.3, 4.3). The received carrier frequency is then effectively 
multiplied by the turnaround ratio G and becomes the 
exciter for the downlink transmitter. For S-band up and 
down frequencies, the ratio G = 240/221, while for S-band 
up and X-band down (S/X) frequencies, G = 880/221. 
The output of the uplink carrier-tracking loop has two 
phase-modulation components: (1) the in-band random 
phase noise, which is a function of the tracking loop 
signal-to-noise ratio p, and (2) components of the uplink 
modulation (e.g., command), which fall within the track- 
ing bandwidth. The intensity of this latter modulation is 
dependent on the uplink modulation index and the modu- 


lated carrier spectrum relative to the tracking loop band- 
width as well as the large frequency characteristics of the 
low-pass loop filter. 

When the phase-modulated output of the uplink track- 
ing loop is frequency-multiplied by the turnaround ratio, 
the phase deviation of the two modulation components is 
also increased by G. The resultant suppression of the 
downlink carrier due to the unwanted modulations, par- 
ticularly for the S/X configuration, can become significant. 
The use of an uplink SCR minimizes to the greatest 
extent possible the turnaround modulation effects. First, 
since p is larger compared to the discrete carrier receiver, 
the random phase noise will be as small as achievable. 
Secondly, one of the virtues of a carrier regenerating loop 
is that it is basically an anti-modulation device in the 
sense that very little of the direct carrier modulation is 
tracked by the loop. This is particularly true of the data- 
aided loop (DAL) introduced in the next section. Thus, 
for all practical purposes, the phase-modulation compo- 
nent due to the uplink sideband frequencies in the vicinity 
of the carrier is eliminated. 

A more detailed discussion of two-way SCRs, along 
with the modulation-induced downlink phase modulation 
components, is presented in Section VII-D. 

C. Known Types of Suppressed Carrier 
Receivers 

The purpose of a suppressed carrier regenerating loop 
is to acquire and track a discrete frequency component 
which is derived in some nonlinear manner from the 
received suppressed carrier signal. The nonlinear opera- 
tion must be of even order and, for the case of single- 
channel signals, is second-order. Three well known regen- 
eration loops will be discussed here: (1) the squaring loop, 
(2) the Costas loop, and (3) the data-aided loop. 

For the purpose of the following discussion, the received 
signal that is input to the receiver will have the form 

x ( t ) = V 2 P D d ( t ) cos ta c t + n, (t) (2-12) 

where the input noise n, ( t ) is bandpass with one-sided 
noise spectral density N 0 W/Hz. (See Eq. 3-2 for a more 
detailed definition.) The suppressed carrier signal com- 
ponent of x ( t ) has the power density spectrum shown in 
Fig. 2-2 

The basic squaring loop is shown m Fig. 2-3. The signal 
x (#) is bandpass-filtered and squared (multiplied by itself) 
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Fig. 2-3. Squaring loop receiver 


to produce z ( t ), which has a discrete frequency compo- 
nent at 2w 0 . The loop locks to this component, and the out- 
put of the voltage-controlled, oscillator (VCO) is divided 
in frequency by two and phase-shifted 90 deg to provide 
the reference r 2 (t) which coherently demodulates the 
carrier. The output ^of the data detector is the maximum- 
likelihood estimate d ( t ) of the transmitted symbols. 

The performance of the squaring loop, as measured by 
the value of L R , is largely a function of the SNR into the 
square-law device (see Ref 2-2) This loss is given by 

L ‘ = [ 1 + t]" (M3) 

where K L is related to the autocorrelation function of 
the input bandpass filter (BPF), H ( s ), and p, is the SNR 
measured in the noise bandwidth of the input BPF. As 
an example of a typical value for L R , suppose that the 
input BPF is a two-pole Butterworth and that the 3-dB 
bandwidth of the filter is set at 4 £R S . Assuming a nomi- 
nal symbol error probability of 1 X 10~ 2 (P D T 3 /N 0 = 2 7), 
p, = 0.608. For the two-pole Butterworth, K L — 0.75, thus, 
L r = 0.62 or —2.1 dB. Even higher losses than this will 
occur if the data stream is coded, since P D T S /N 0 <<27 


We see that the squaring loss is a function of input SNR, 
which in turn depends upon the bandwidth of the input 
filter. Obviously, in order to minimize the loss, the filter 
bandwidth should be made as narrow as possible without 
attenuating a significant portion of the signal power. This 
presents certain mechanization problems. First, the center 
frequency of the BPF may be relatively high (generally 
the last IF frequency in the receiver), while the symbol 
rate £R e may be quite low, requiring that the filter be very 
high Q. This not only introduces linearity and stability 
problems but can also result in high costs. Furthermore, 
each time the data rate is changed, it is necessary to 
change the filter, further increasing costs, and introducing 
a degree of inflexibility. Thus, the use of the squaring loop 
for suppressed carrier, multiple-data-rate applications has 
severe limitations It has been used for the tracking and 
demodulation of low-data-rate, biphase-modulated sub- 
carriers of the type used in the Mariner 1969 (and follow- 
ing) two-channel command systems (Ref. 2-3). 

The Costas loop, shown in Fig 2-4, does not require 
that the input BPF have a specified bandwidth, rather, its 
performance depends upon the specification of the two 
low-pass filters, G (s). Actually, the Costas loop is nothing 
more than a low-pass version of the squaring loop (see 
Ref 2-4 for bandpass-low-pass equivalents), and has the 
identical stochastic differential equation as that of the 
squaring loop, provided that the noise spectrum formed 
by the filters for each configuration is the same (Ref. 2-2) 
Since it is generally easier to build switchable, active, low- 
pass filters, the Costas loop has a mechanization advantage 
over the squaring loop for multiple-data-rate applications 



Fig. 2-4. Costas loop receiver 
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It should be noted that for both the squaring and Costas 
loops, no attempt is made to optimize the receiver per- 
formance with respect to the modulating signal d (t). That 
is, the nature of d (t), other than its power spectrum, plays 
no real part in the formation of the loop error signal. An 
optimum receiver would, however, make proper use of the 
data estimate d (t). A receiver configuration which accom- 
plishes this is the data-aided loop, the basic configuration 
is shown in Fig 2-5 The topology of the DAL is very 
similar to that of the Costas loop, there are, however, 
major differences. First, the data detector is now an inte- 
gral part of the loop. Secondly, a time delay equal to the 
symbol duration T s = 1/£R, is employed. Thirdly, the low- 
pass filters G (s) are not required, and the input BPF serves 
only to limit the total input noise voltage. 

The performance analysis of the DAL is the main sub- 
ject of this report, beginning with Section III. Drawing 
from one of the results presented in Section III-C, we may 
compare the DAL performance to that of the squaring 
(or Costas) loop. For the case where the symbol rate fR, 
is large compared to the loop bandwidth W L , the efficiency 
loss for the DAL is given by 


Finally, it should be pointed out that the DAL offers 
certain mechanization advantages over, say, the Costas 
loop. It is very amenable to simple all-digital implementa- 
tion, which makes the problems of data-rate switching 
virtually disappear. Details on mechanization are given in 
Section IV. 

The data-aided loop as described in this report has yet 
to find application in any operational system. There are, 
however, two closely related mechanizations which are 
presently being employed. 

The first application is found in the all-digital command 
detector (Ref. 2-6), where the subcarrier tracking loop is 
indeed a form of the DAL The second is the subcarrier 
tracking loop found in the DSN’s subcarrier demodulator 
assembly (SDA) (Refs 2-7 and 2-8) In the SDA, the loop 
is actually closer to a Costas loop than a DAL, with the 
exception being that a hard limiter follows the low-pass 
filter (LPF) which contains the demodulated data sym- 
bols, thus making a continuous, nonintegrated estimate of 
the data symbol polarity 


L* = ( 1 - 2 P‘) 


(2-14) 


where P’ e is the average symbol probability of error. Using 
the previously cited conditions for the squaring loop, we 
readily find that L R = 0.98 or —0.1 dB. Thus, for the stated 
example, with respect to loop signal-to-noise ratio, the 
DAL performance exceeds that of the squaring loop by 
2.0 dB. A more parametric comparison of these loops is 
given in Ref. 2-5. 



Fig. 2-5. Suppressed carrier data-aided receiver 


D. Differences Between DAL Suppressed 
Carrier and Discrete Carrier Systems 

For discrete carrier systems, the RF tracking loop need 
only acquire and track the carrier line in the received 
spectrum. The tracking loop has no knowledge of the data 
modulation The carrier loop does, of course, provide ref- 
erences for demodulating the data-modulated subcarrier 
from the carrier. 

The DAL, on the other hand, must have credible esti- 
mates of the data symbols before a useful tracking error 
signal can be formed. However, the recovery of data 
symbols implies coherent earner demodulation. Thus, the 
acquisition of the DAL would seem to involve one of those 
incredible schemes to determine “which comes first, the 
chicken or the egg.” Actually, this is not the case, and it is 
quite easy to obtain good data symbol estimates even 
though the loop is very much out of lock due to the fact 
that the symbol transitions are not affected by carrier 
lock, thus permitting the data detector to synchronize and 
function The effect of the carrier being out of lock is to 
amplitude-modulate the data symbol stream by the loop 
error signal 

All suppressed carrier regeneration tracking loops, in- 
cluding the DAL, have double-frequency loop error sig- 
nals with respect to those of the discrete carrier loops 
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This gives rise to differences in the stochastic performance, 
and introduces a data phase ambiguity. 

The DAL also has a finite delay function which can con- 
tribute to instability under certain conditions. In terms of 
attendant functions, the DAL has an envelope type of 


automatic gain control (AGC), and an absolute value type 
of lock detector (see Section V). 

The differences between the DAL and discrete tracking 
loops are summarized in Table 2-2, along with the text sec- 
tions in which detailed discussions are found. 


Table 2-2. Differences between discrete carrier and data-aided tracking receivers 


Item 

Discrete earner loop 

DAL 

Sections 

Loop operation dependent on data modulation 

No 

Yes 

II-D, m-A, III-B, III-C, IV-B 

Subcarrier required for single-channel 
modulabon 

Yes 

No 

II-A 

Form of linear loop error signal 

0 

20 

III-C 

Data estimate has phase ambiguity 

No 

Yes 

m-c, vi - a 

Delay function instability 

No 

Yes 

III-D 

AGC system 

Coherent amplitude 

Envelope 

VI-A 

Magnitude of modulation feedthrough 

Spectrum, loop bandwidth, and 
modulabon mdex dependent 

Essentially zero 

II-B, VII-G 
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III. A General Theory for the Data-Aided Receiver 


A. The Basic Concept and Its Functional 
Configuration 

The basic idea of a data-aided receiver (DAR) centers 
around using the power in the composite received signal 
(carrier and data-modulated subcarriers) sidebands to 
enhance the signal-to-noise ratio in the bandwidth of the 
carrier-tracking loop The vehicle through which this is 
accomplished employs the principle of decision-directed 
feedback, and the carrier-tracking loop to which it is 
applied is referred to as a data-aided loop (DAL). The 
underlying idea behind the data-aided approach is that 
an enhancement of loop signal-to-noise ratio (SNR) not 
only improves the carrier-tracking function (by reducing 
the loop’s mean-squared phase noise and average rate of 
cycle slipping) but also results in an improvement in the 
data detector performance because of the accompanying 
reduction in noisy reference loss. 

For a single-channel system wherein a discrete carrier 
is phase modulated by a biphase-modulated data sub- 
carrier, the sideband power being centered around the 
subcarrier frequency would ordinarily be filtered out by 
the carrier tracking loop [eg.,a phase-locked loop (PLL)], 
and hence does not affect carrier-tracking performance. In 
the DAR, an estimate of the biphase-modulated data sub- 
carrier is formed which, when fed back to the DAL, can 
be used to recover the power in the sideband components 
for carrier-tracking purposes 

More specifically, consider the block diagram illustrated 
in Fig. 3-1. The upper half of the configuration constitutes 
a PLL used in present-day tracking receivers. The modi- 



Fig. 3-1. Data-aided carrier-tracking loop 


fication proposed here makes use of a quadrature channel 
to insert a signal component into the voltage-controlled 
oscillator (VCO), whose power is proportional to that 
of the biphase-modulated data subcarrier. One principal 
point to note is that the shape of the equivalent phase- 
detector characteristic (S-curve) is affected by the power 
m the signal’s sidebands, the phase jitter in the subcarrier 
tracking loop, the symbol sync jitter, and the conditional 
probability of error of the data detector. However, cursory 
examination of these effects indicates that the symbol sync 
jitter and subcarrier sync jitter are small relative to the 
radio loss caused by the carrier-tracking loop phase error; 
consequently, essentially all of the sideband power can be 
recovered and used for improving the principal source of 
system degradation, viz., the noisy reference loss. 

In Section II, we established a case for suppressed 
carrier signaling, that is, where the data are directly 
biphase-modulated on the carrier, thus, we shall focus our 
entire attention on the specific version of the DAL re- 
ferred to as the suppressed carrier DAL (denoted hence- 
forth simply by DAL for the sake of brevity) (see Fig. 3-2). 

B. Development of the Stochastic Integro- 
Differential Equation of Operation 

In this section, we develop the loop equation of opera- 
tion to the extent that it applies to Fig. 3-2. For a single- 
channel suppressed carrier system, the received signal. 



Fig. 3-2. Suppressed carrier data-aided carrier-tracking loop 
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x ( t ), is of the form 


where the signum function sgn x is defined by 


x ( t ) = V 2 P D sin + yrf(t) + 0J + n,(<) 

— V 2 P D d ( t ) cos ( t ) + n, ( t ) (3-1) 

where 6 = Q 0 t + 0 o is a slowly varying function of time 
which is related to the doppler offset Q 0 , the function $ (t) 
is defined by $ (#) = a c t + 0, and the random process 
(n, (t)} is assumed to be narrowband and gaussian. A sam- 
ple function representation of this process is given by 

n, ( t ) = V2 [n, ( t ) cos # ( t ) — n 2 ( t ) sin $ (?)] (3-2) 


(+1 x > 0 

j-1 x<0 


(3-7) 


Expanding x (<) of Eq. (3-1) into its signal and noise com- 
ponents, we can rewrite (3-6) as 


et(t ~ T,) = sgn |K,K lm VPd f T ‘d (t) cos <f> (t) dt 


+ 


(V [t,4>(t)]dt | 


(3-8) 


where n, (t) and n 2 (t) are independent gaussian random 
^processes with single-sided spectral densities N 0 , W/Hz 
(It is convenient to characterize the reference signals r,(t) 
and f 2 (t) in Fig. 3-2 by 


fi ( t ) = V~2 K, cos'S' (t) 
r 2 (t) = -Y2K 2 sin^(t) 


(3-3) 


where '3> (t) — <a c t + 'O', and 'O' is the loop estimate of 0 
Assuming a lower arm multiplier gain of K 2m and neglect- 
ing double-frequency components, the output z ( t ) of the 
filter F (p) is given by 


z(t) = K 2 K 2m F(p) exp [ — pT„] 

X (V Pod ( t ) d (t) sin <j> (t) + d ( t ) N 2 [t, <£ (t)] } 

(3-4) 


where 


A 

IV, [t, <f> (t)] — n, (t) cos 4> ( t ) — n 2 (t) sin (t) (3-9) 

is also a baseband white gaussian noise process of single- 
sided spectral density 1V 0 W/Hz. Furthermore, IV, [t, <f> (t)] 
and N 2 [t, <f> (#)] are uncorrelated Also, for simplicity of 
notation, we omit hereafter the dependence on t in all 
variables which are functions of time 

When the two-sided loop bandwidth W L of the DAL is 
small relative to the data rate £R, = 1 /T,, then the phase 
error <p can be considered to be constant for many T« 
second intervals Under these conditions, the conditional 
error probability of the data detector, i.e., the probability 
that d^d, given a phase error of </> rad, is essentially 
independent of t and given by 


where 


P ’’ (<t>) = jr erfc m cos <f>) (3-10) 


N 2 [t, <j> (f)] = n, (t) sin <j> (t) + n 2 (t) cos <t> ( t ) (3-5) 


where 


is a baseband “white gaussian noise process” of single- 
sided spectral density N 0 W/Hz, and 


R s = 


PoT, 

N„ 


(3-11) 


= *(*) -$(t) 

\ 

is the loop phase error. 

\ /\ 

The upper arm of Fig. 3-2 forms the data estimate d 

Assuming a'multiplier gain of K im for the upper arm cross- 
correlator, then 

^(t - T„) — sgn ’x(t)r(t)dt~^ (3-6) 


and 


erfc: 


A _2_ 

“vrj. 


exp ( — t 2 ) dt 


(3-12) 


Furthermore, insofar as it^ effect on the loop equation is 
concerned, the product dd can be replaced by its statisti- 
cal average, i.e., E {dd} = ( + 1) X Pr (d =ct) + ( — 1) X 
Pf {d ^=a) — { + 1) X [ 1 — P® (c£)] + ( — 1) X Pj (<£) = 1 — 
2P® ( 0 ). The reason for this is that the loop filter being 
low-pass in nature will tend to extract the time average 
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(or dc value) of the random signal ad, and since dd is 
ergodic, its time and statistical averages are equal Using 
these assumptions, Eq. (3-4) reduces to 

z = K 2 K 2m F ( p ) exp [- pT s ] 

X ("V P D erf (VTC cos <£) sin <f> + dN 2 (f, 0)} 

(3-13) 

A 

where erf x = 1 — erfc x. 

It is clear that, in general, the transfer function factor 
exp (— pT s ), with p = affects loop stability and reduces 
the signal acquisition or pull-in range as discussed in 
Sections D and E. For the moment, we shall ignore these 
problems and make a simplifying assumption which en- 
ables us to neglect this exponential factor in regard to 
predicting steady-state performance If W L T a < < 1 (i e , 
the usual case of interest), then exp ( — jo>T a ) is approxi- 
mately unity for all <o within the bandwidth of the DAL. 
Hence, from a steady-state performance standpoint, this 
factor has negligible effect Mechanistically, however, the 
delay in the lower Joop is important in assuring that d is 
multiplied by the d corresponding to the same symbol 
interval. 

Making the usual assumption that the instantaneous 
VCO phase 0 is related to its input z through 


of the results which follow to the second-order DAL is 
straightforward based upon the contribution made in 
Ref. 3-1 The essential details of this extension will be 
presented at the end of our discussion 

The steady-state performance of any carrier-tracking 
loop operating in its nonlinear region can be characterized 
by (1) its phase noise o% (i.e., the variance of the modulo-27r 
reduced phase error) and (2) its average rate of cycle slip- 
ping S. With respect to calculation of the first of the above 
two performance parameters, it is essential that we char- 
acterize the probability density function (p.d.f.) of the 
modulo-27r reduced phase error, which for simplicity we 
shall also denote by <f> Hereafter, when we use the term 
phase error, we shall implicitly mean the modulo-27r re- 
duced version unless otherwise specified 

1. Steady-State Phase Error Probability Density Function 

Using results given in Ref. 3-1 (Eq. 32), it is easy to 
show that the p.d f of the phase error process is given by 

p(4>) = N exp [~U 0 (<},)] / exp [U 0 (x)] dx, 

h 

(3-16) 

where N is the normalization constant and the potential 
function U 0 (<f>) is given by 



(3-14) 


where K v is the gam of the VCO, we find, from Eqs. (3-13) 
and (3-14), and the fact that <£ = 6 — 6, 


U 0 (<t>)= — o<t> + KV Pd [cos </> erf (V R s cos </>) 

K-oo ( 

+ y^ exp(-R t cos 2 <ft)]| (3-17) 


<j> = 6 — KF (p) [ V Pd erf (V R* cos <j>) sin rf> + dN 2 (t,<}>)] 

(3-15) 


with 




VqK 2 

2 


(3-18) 


where K — K 2 K 2m K v represents the loop gain Equa- 
tion (3-15) represents the stochastic mtegro-differential 
equation of loop operation for the configuration illustrated 
in Fig. 3-2. 


The loop gain K can be related to the system parameters by 


K = 


4 

sVP^erfVTL 


(3-19) 


C. Steady-State Performance 

It is simplest at the outset to consider the steady-state 
performance characteristics of the first-order DAL, where 
the loop filter F (p) is of the form F (p) = 1 The extension 


8 


a 2 

~ w l t„ 


(3-20) 
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A VP^KerfVR* 


(3-21) 


is the two-sided bandwidth of the DAL. Substitut- 
ing (3-18), (3-19), and (3-20) into (3-17) gives, upon 
simplification, 


Vo ( <t > ) = - ~ R ° s erf YR~, 


j^cos^erf (V R s 


cos <j>) 


+ y exp ( - R s cos 2 <j>) J (3-22) 


In terms of the loop signal-to-noise ratio, 


A 4 V~Rd erf VRl 
p N 0 K 


_ 2P D erf 2 V~Rr 
N„W t 

= R s S erf 2 V R s 


(3-23) 


and the steady-state phase error (often called static phase 
error or SPE), 


+'* _ Sin_1 [ YP^K erf Yk] (3 ' 24) 


Eq. (3-22) may be rewritten as 


U 0 (<j>) — — (p sin <j> se ) <j> — p |” 


cos $ erf (V Rs cos <£) + (1/V ttR 8 ) exp ( — R s cos 2 $) ' 
erf V R, 


We note that when n 0 = 0, Eq. (3-25) reduces to 


U 0 (<t>) — ~p 


_ cos (/> erf (YK cos <j>) + (1/V ttR,) exp ( — R, cos 2 <ft) 


erf \TiL 


(3-26) 


For a first-order PLL, the potential function U n (<f>) is simply 


U o (<£) = — p COS <t> 


— p j s\n <f> d<f> 


(3-27) 


Thus, the effective loop nonlinearity (or S-curve) of the DAL is given by 


d_~ 

d<j> _ 


d f cos <j> erf (V R„ cos <j>) + (1/V 7rR s) exp ( — R, cos 2 <f >) " 


erf V 8? 


of (tf>) — 


erf (V R s cos <£) sin<£ 

erf V~R» 


(3-28) 


and is plotted vs <j> in Fig. 3-3, with R, as a parameter. We 
note that in the limit of infinite data signal-to-noise ratio, 
(3-28) becomes 


where 5, (<j>) is a unit power square wave in <f> with period 
equal to 2 n, i e., 


5,(*) = 


+ 1 — 2 — 1“ 2n?r — <f> — — f 2n- 


- 1 — — h 2n— — 0 — -g — I - 2n7; 


of (4>) = sgn (cos 4>)sin<f, — S q (<j>) sin </> (3-29) 


n — 0, ±1, ±2, 


(3-30) 
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nize the phase of unique sync words, or (3) for uplinks 
(ground to spacecraft), use the procedure outlined in Sec- 
tion VI- A. 


Figure 3-4 illustrates the phase error pdf., for 

Q 0 = 0, 8 = 10, and R e as a parameter. Notice from this 
figure that the DAL exhibits the same 180-deg phase am- 
biguity as the squaring or Costas loop, because the stable 
lock points are separated by only 180 deg. This ambiguity 
can be resolved in several different ways: (1) differentially 
encode and decode the symbols, (2) search for and recog- 


2. Phase Noise Performance 

Assuming that the phase ambiguity is perfectly resolved, 
then the phase noise <4 can be computed numerically from 

/ V2 r /•*/ 2 -|2 

#>(<#>) <4 J (3-31) 


For R, — » 0, we get 




Fig. 3-3. Equivalent loop nonlinearity (S-curve) of the DAL 
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Fig. 3-4. Probability density function of the loop phase error 


where p ( <j > ) is given by (3-16) together with (3-25). When 
R, becomes large, a\ behaves as follows: 


4 J 

7T/2 

4> 2 exp (p cos <j>) dd> 

0 


rir/2 

/ exp (p cos 4>) d$ 

0 


which, for large p, becomes 


(3-32) 


In comparing the carrier-tracking efficiency of discrete 
and suppressed carrier systems, it is important that the 
comparison be made at a fixed value of total transmitter 
energy-to-noise ratio P t T s /N 0 and symbol time-bandwidth 
product 8 For the suppressed carrier system (employing 
a DAL for carrier tracking), the data signal- to-noise ratio 
R« is identical to P t T s /N 0 . Thus, the loop signal-to-noise 
ratio p is a direct function of 8 and P t T s /N 0 through 
Eq. (3-23). For the discrete carrier system (employing a 
PLL for carrier tracking), the loop signal-to-noise ratio 
p — 2 P c /N 0 W l is related to P t T s /N 0 , R«, and 8 by 



(3-33) 

Also, 


Thus, insofar as phase noise is concerned, the first-order 
DAL has identical behavior in its linear region of opera- 
tion as the PLL, This result is in contrast to the other types 
of suppressed carrier-tracking loops discussed in Sec- 
tion II, whose phase noise in the linear region is degraded 
by the squaring loss. 


P = [(PtT,/N o) - R S J 8 (3-34) 




-r 


, , exp (p cos <f>) 

* HW # 


00 



k = l 


(3-35) 
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Figures 3-5 and 3-6 plot <4 vs P T T a /N 0 for fixed values 
of S = 10 and 8 = 50, respectively. For the discrete carrier 
system, o% is computed from Eq. (3-35) together with 
(3-34), where R a has been set equal to unity and the range 
of P t T s /N o is selected to correspond to a range of values 
for p of 0 to 10 dB. The value of R„ was selected to corre- 
spond to the minimum acceptable error probability per- 
formance in the presence of a perfect carrier reference 
For the suppressed carrier system, o% is computed from 
Eqs. (3-31) and (3-16) and plotted over the same range of 
P r T e /N 0 as the curves for the discrete carrier case We 
note from the figures that at each value of P T T a /N 0 , the 
DAL has a smaller <4 than the PLL. Furthermore, since for 
the PLL, the relation for R s given by (3-11) implies a zero 
value of loop signal-to-noise ratio (see Eq. 3-34) or equiva- 
lently, an infinite value of <4, then for any finite value of <j% 
and a fixed value of P t T s /N 0 , the suppressed carrier sys- 
tem allows for a larger value of R„ than the discrete carrier 
system. This increase in data signal-to-noise ratio mani- 
fests itself in an improved error probability performance. 



Fig. 3-5. Variance of the loop phase error as a function of 
total energy-to-noise ratio (6 = 10) 


In conclusion, the suppressed carrier system with a 
DAL has a better error probability performance than the 
discrete carrier system with a PLL at any given value of 
transmitter energy-to-noise ratio because (1) the carrier 
reference has less jitter, resulting in a smaller noisy refer- 
ence loss, and (2) the data signal-to-noise ratio is larger 

3. Average Probability of Error of the Data Detector 

When the loop phase error can be considered essen- 
tially constant over several symbol intervals, then the con- 
ditional probability of error of the data detector is given 
by Eq. (3-10). The average error probability (assuming 
that the 180-deg phase ambiguity is perfectly resolved) 
is then obtained by averaging P* (<f>) over the p d.f. p (<j>) 
of the phase error as given by (3-16), i.e., 


/ ir/2 

P 8 e (</>) P (<t>) <U> 

ir/2 



Fig. 3-6. Variance of the loop phase error as a function of 
total energy-to-noise ratio (S = 50) 
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4. Average Rate of Cycle Slipping 

A general expression for the average rate of cycle slipping S of a first-order suppressed carrier-tracking loop is (Ref 3-2) 


— 1 , 7TQo 

S ~ — CoDoKoo coth — 

u A-oo 


(3-36) 


where K 00 is defined by (3-18), and 

c - - [/ 


A C rir/a / fP 

C 0 = / exp [ — U 0 (<£)] < 1 + D 0 / exp [U 0 (x)] dx 

^J-V/2 V J-TT/2 


KT 


A exp [-£/<,( -tt/ 2)] - exp [-Up (tt/2)] 

r v/i 

exp [-C7 0 (tt/2)] / exp [f7 0 (<f>)] d<f> 

J -V/ 2 

Since the potential function U 0 (</>) can always be expressed in the form 

On f 

Uo (4>) = - <t> + p JJ(4>)d+ 

then, after much simplification, 


(3-37) 


(3-38) 


K.q o cosh {ttCIo/Koo) 


r ir/2 r ir/2 r Tt/2 r <t> 

exp ( ttOo/Koo) / exp [-U 0 (<t>)]d<t> / exp [t/ 0 (4>)] d<j> ~ 2 sinh (irflo/Koo) / exp [-U 0 (<j>)] dj> exp [17 0 (x)] dx 

J-V/2 J -Tt/2 J -Tt/2 J-tr/2 

(3-39) 


For zero detuning, Eq. (3-39) reduces to 

S = - 


r Tt/2 r 1T/2 

/ exp [ — 1/ 0 (</>)] d<f> / exp [C/ 0 d<t> 

I -Tt/2 J -TT/2 


(3-40) 


Applying (3-39) to the DAL, we get 


s 


(4/p) cosh [(tt/ 2) p sin <j>„] 



w, ~ 

exp [(tt/2) p sin <£„] 

[*/ 2 

f exp [ — U 0 (<£)] d<j> j 

-TT/2 J 

fV/2 

1 exp [U„ (</>)] d</> - 2 swh[(vr/2)p sin $„] i 

-IT/2 J 

f exp [ C7 0 (<#>)] I 

-TT/2 J 

f exp [C/ 0 (x)] dx 

-V/i 


(3-41) 


which, for fi 0 = 0, reduces to 


pj exp £ — P J d (x) dx J d$ j exp p j d (x) dx J d<j> 

= | pj*' 2 exp £ - P J * d (*) dx ] d<j> j V2 exp [ p j' * d (x) dx J dj> j '* 
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Fig- 3-7. Normalized average rate of cycle slipping as a 
function of total energy-to-noise ratio (8 = 10) 


Fig. 3-8. Normalized average rate of cycle slipping as a 
function of total energy-to-noise ratio (8 = 25) 


Figures 3-7 and 3-8 illustrate the variation of S/W L , with 
P t T,/N 0 at fixed 8 for the same parameters as in Figs, 3-5 
and 3-6, respectively. The appropriate equation for the 
PLL is 

WT = ttVo(p) (3 ' 43) 

Also included in Fig 3-7 are some experimental points 
taken from a second-order DAL for the cases of r = 10 
and t = 2. These data are included to verify the slope of 
the S/W L curve for the DAL, as well as to provide an 
assessment of the validity of Eq. (3-36), when applied to 
second-order loops. Noting that the first-order DAL per- 
formance is approached as r-» oo, the experimental data 
(r = 10) strongly confirm the analytical results. 

Asymptotic (large p) approximations to (3-42) and (3-43) 
are possible For the PLL, 

= ~ exp ( - 2p) (3-44) 


while for the DAL, 

w L = exp [ ~ ( 2 ~ t) p ] (3 ' 45) 

5. Extension to Case of Second-Order Loops 

For the second-order DAL with loop filter 

F ( s ) = TTs^7 M 

the p.d.f. of the phase error is still characterized by (3-16); 
however, the potential function of (3-25) must be modified 
as follows: 

UoM = —p<i> 

[cos <f> erf (V~R7 cos <f>) + (1/V ttR s ) exp ( — R, cos 2 </>)] 

_a L irfVRl J 

(3*47) 
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where 


\ T ) rG 2 (<j>) 



X [n 0 -V^KerfVX(l-Fi)GrW] 
+ a ef (</>) 


F, = ^. 

T i 

A f+1 

W, =^r— 


r = VT^Kr^erfVX 


G(<j>) — ef (<]>) — of (< f> ) 


(3-48) 


and the overbar denotes statistical expectation with re- 
spect to the p.d.f. p (</>) We note that the above reduces 
to the results previously given for the first-order DAL by 
letting r approach infinity and F t = 1. 


Using the p d.f. as characterized by the above param- 
eters, the phase noise performance of the second-order 
DAL is still characterized by Eq. (3-31), and its asymp- 
totic performance for large loop signal-to-noise ratio is 
given by (3-33), with p replaced by a. Insofar as the aver- 
age rate of cycle slipping is concerned, Eq. (3-36) can be 
used to approximate the performance of the second-order 
DAL for large r provided that the error, as indicated by 
the experimental data in Fig 3-7 is tolerable. For small r, 
Eq. (3-36) is not very useful, as shown by the experimental 
data for r = 2 in Fig. 3-7 A similar variation with r has 
been shown for the standard PLL in Ref. 3-2 (Fig 10-22) 


D. Effects of the Delay Function 

Unlike the PLL or the other types of suppressed carrier 
loops discussed in Section II, the DAL requires an ideal 
delay element (see Fig. 3-2 and Section III-B) in its lower 
arm for proper operation. The primary purpose of this 
delay element is to delay the signal emanating from the 
lower arm cross-correlator by an amount equal to one 
symbol time T g so that in the multiplier which follows, 
a particular data symbol is multiplied by the estimate of 
that same symbol; i.e., the matched filter in effect intro- 
duces a delay of one symbol time Aside from its neces- 
sary value, the presence of a delay element in a closed- 


loop configuration produces several deleterious effects on 
the loop’s performance. The first of these to be considered 
is the effect of the delay element on loop stability. 

1. Effect of Delay Element on Loop Stability 

In the past, many techniques (e.g., the Nyquist crite- 
rion, root-locus diagram, Routh-Hurwitz criterion, phase- 
plane, Lyaponov function) have been developed for 
performing a stability analysis on a closed-loop system 
Unfortunately, all of these available techniques are appli- 
cable only in a noise-free situation In the case of the 
second-order PLL, this restriction on applicability does 
not present a problem since a linear, noise-free analysis 
predicts unconditional stability, whereas the nonlinear 
techniques demonstrate the existence of an adequate 
pull-in range For the DAL, the problem of loop stability 
becomes more poignant because of the presence of a delay 
element in the open-loop transfer function. This delay, 
which is ideally equal to the reciprocal of the data rate, 
will, at sufficiently low data rates, cause instability prob- 
lems even on a linear, noise-free basis. Furthermore, the 
nonlinear stability problem is complicated by the fact that 
the loop S-curve of (<f>) is now a function of signal-to-noise 
ratio, and hence a noise-free analysis does not tell the 
whole story To get around this problem, we shall con- 
sider only the nonlinear stability problem in the absence 
of noise, in which case the loop S-curve is given by 
Eq (3-29). The results obtained from this special case can 
then be thought of as the “best possible” nonlinear sta- 
bility performance and, together with those obtained from 
the linear analysis, shall be used to provide an approxi- 
mate answer to the question of how low in data rate one 
can go before the DAL becomes unstable. 

The principal tool used in the linear stability analysis 
which follows is the Nyquist diagram (and the accom- 
panying Nyquist criterion) as applied to the open-loop, 
noise-free transfer function of a generalized second-order 
tracking loop. For the nonlinear stability performance, we 
resort to a method based upon an approximation to the 
phase-plane behavior 

a. Linear Stability Analysis. Consider a generalized 
phase-locked tracking system whose noise-free differential 
equation of operation is given by 

* = (3-49, 

where 8 is the input phase to be tracked, <f> = 6 — 8 is the 
loop phase error, is the loop’s estimate of 8, F ( p ) is 
the loop filter, of ( <j> ) is the system nonlinearity normalized 


18 


JPL TECHNICAL REPORT 32-1583 



to unit slope at the origin, AK is the open-loop gain, and 
T d is an arbitrary ideal delay. For the DAL, ef (<f>) is 


defined in (3-28), A = \TfW \[R S> and T d = T e , i.e, the 
reciprocal of the data symbol rate. The open-loop transfer 
function T ( p ) = 'd'(p)/& (p) is from Eq. (3-49). 


T(p) = AK 


EM 

p 


exp ( — pT s ) 


(3-50) 


then the second-order tracking loop will be absolutely 
stable m the linear sense. Assuming fr, > > r 3 , the solu- 
tion of Eq. (3-55b) for x 0 is 


*0 > 



1 + V 1 + 4/r 2 
2 


(3-56) 


and hence, 


where ef (</>) has been linearized to 4> Substituting (3-46) 
into (3-50) and letting p = jo>, we get 


tan* 1 x n 


o 


(3-57) 


m,. S AK( 1+7o.T 2 ) , . r„s 

r(; “ ) = 7 < ;(iT/ <!> rT) exp(_;wTs) ( 3 - 51 ) 

A 

Introducing the normalization x = <ut, and letting r = 
AKtI/t i, Eq. (3-51) simplifies to 


T (/*) = 


rQ-i/rg) 2 [1 + /x^/rQ] exp [ — ;3c(T,/n)] 
jx(l + jx) 


(3-52) 


The quantity r is related to the loop damping £ at zero 
delay by 


r = 4£ 2 


(3-53) 


Substituting expressions (3-56) and (3-57) in Eq (3-55a) 
gives 


T s , P.V (tan- 1 [f V (1 + V 1 + 4/rQ/ 2]} 
T2 r V (1 + V 1 + 4/r 2 )/2 


(3-58) 


where P.V. denotes the principal value. Equation (3-58) 
then represents a bound on normalized loop delay in order 
to guarantee linear stability. As an example, for r = 2 (i e , 
£ = 0.707), T,/ t 2 must be less than 0.52 (see Fig. 3-9) 


It is desired to study the magnitude and phase charac- 
teristics of T ( pc ) as defined in Eq. (3-52) and from these, 
using the Nyquist criterion, determine the set of normal- 
ized frequencies { x „ ) at which the phase characteristic, 
arg T ( jx ), equals (2 n — 1) v rad, n = 0, 1, 2, . Thus, 

the set {*„} is the solution to 


The quantitative results given above rely heavily on the 
ability to make the small phase error assumption which 
allowed us to use the linear model and a linear stability 
criterion. In fact, we notice that the linear stability analy- 
sis determines only whether the steady-state phase error 


tan- 1 £ x n J ~ tan_1 - x n ~ = (2 n - 1) tt 

(3-54) 

In particular, the solution of Eq. (3-54) for n = 0 (i.e , the 
first crossing of the horizontal axis by the Nyquist plot) 
must have a corresponding magnitude of T (jx) which is 
less than unity in order to ensure stability. Thus, if the 
delay T s is such that 


— tan’ 1 Xq + tan- 1 x 0 + x 0 — = y (3- 


and 


55a) 


»*(*7t,H1 + *8(t,/t 1 )*] 

= 4TTT4) <' < 3 ' 55b) 



Fig. 3-9. Nyquist piot for second-order tracking loop with delay 
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becomes zero or unbounded and hence does not predict 
the existence of stable and unstable limit cycles Thus, one 
should regard the above results as a qualitative guide 
(possibly even as a loose upper bound) to the real behavior 
which must be predicted by a nonlinear model. 


b. Nonlinear Stability Analysis. Despite the fact that we 
must simplify the nonlinear stability problem by assuming 
a noise-free situation— i e , the loop S-curve is given by 
(3-29)— an analytical solution for the critical value of delay 
above which the loop becomes unstable is difficult to come 
by In fact, the only solutions which have thus far been 
obtained for problems of this nature pertain only to the 
PLL. More specifically, for a PLL with a loop filter given 
by (3-46) and an arbitrary delay T d in the loop, the abso- 
lute instability region is given by (Ref. 3-2, Eq 10-67) 


T d 

/ 19.34t 2 a 

> A 

/ 2 7 i C 

T 2 ^ 

Tjf 2 


(3-59) 


In other words, for values of T d < t 2 ^ c , the loop may or 
may not reach a stable lock point, depending, of course, 
on the initial frequency offset. If the initial frequency 
offset is too large, then the loop goes into a stable limit 
cycle These considerations on the loop’s behavior in its 
region of conditional stability will be discussed in Sec- 
tion E. 

It is interesting to note that the region of absolute insta- 
bility as described by (3-59) also applies to a tracking loop 
with a sin 2 <f> characteristic, e.g , a squaring loop Thus, 
despite the fact that we cannot find an expression for the 
absolute instability region of the noise-free DAL, we 
realize that the S-curve of (3-29) is also a nonlinear func- 
tion of 2$, and on this basis alone, we might conjecture 
that (3-59) applies at least in a qualitative sense Further 
results of a quantative nature must at this point come from 
experimental evidence. Experimentally it has been found 
that for a DAL with r = 2 and ti/t 2 = 200, a stable situa- 
tion is achieved when > 2 W L , where W L is the loop 
bandwidth at zero delay as given by (3-21). Since (R, — 
I /T a and T d =T,, then expressing the inequalities of (3-58) 
and (3-59) in terms of £R S and W L , the corresponding in- 
equalities necessary to guarantee stability become 


> l.3W t (linear theory) 
cR a > 4.3 W L (nonlinear theory) 


(3-60) 


2. Effect of Delay Function on Loop Bandwidth and 
Natural Frequency 

In order to assess the effect of the delay function on 
loop bandwidth and natural frequency, it is essential that 
we characterize the linear model of the DAL in terms of 
its closed-loop transfer function Letting of (</>) be approxi- 
mated by <j>, then the closed-loop transfer function which 
relates the phase estimate 6 to the input phase 8 is 


AKF (s) exp (— sT 8 ) 

5 9 (s) s + AKF ( s ) exp ( — sT „ ) 


(3-61) 


where we have set A = V Pd erf \TRs. Substituting (3-46) 
for F (p) into the above yields, upon simplification. 


H(s) 


1 + S To 


1 + {r 2 + [exp (sT s )]/AK) + {[r,exp (sT,)]/AK} s 2 

(3-62) 

Letting 

S St 2 

V = T k /t 2 (3-63) 


then using the definition of r as given in (3-48), we have 
the normalized form 


H (s') = 

1 + s' 

1 + {1 + [exp (s'r)]/(fT, / r 2 )} + {[exp (, sT')]/r } s' 2 

(3-64) 


If rn > > t 2 (the usual case of interest), then 


H(s') 


1 +s' 

1 + s' + {[exp (sT')]/r} s' 2 


(3-65) 


The linear loop bandwidth W L is defined by 

W L = f ’ X H (s) H (— s) ds (3-66) 

or equivalently, 

W L r 2 = -J— f ’ X H (s') H (—s') ds' (3-67) 

" 7r / J -JX> 

Letting s' — /V, 

w l t 2 =~J X imh'WcL' (3-68) 
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Substituting (3-65) into (3-68), we get 



1 + jta 

[1 - (c/ 2 /r) cos u'T'] +/</[! + (to' sin o/T')/r] 


diJ 


or 


!/■“ j / 

t J 0 1 + */ 2 [1 - (2/r) cos o,T'] - <„' 3 [(2/f) sin </T'] + o/yr 2 


(3-69) 


At zero delay, from (3-48), 

W, T2 = (3-70) 

Plotted in Fig 3-10 is the ratio of W l t 2 as given by 
(3-69) to its zero-delay value (Eq. 3-70) vs normalized 
delay T' = T„/t 2 , with r as a parameter We note from the 
figure that the presence of delay in the loop can signifi- 
cantly increase the effective linear loop bandwidth when 
the value of this delay (reciprocal of the data rate) relative 
to the filter time constant r 2 becomes appreciable. 

The natural frequency of the linear model of the DAL 
(or for that matter, any carrier-tracking loop) can be mea- 
sured as follows (Fig. 3-11). A carrier (at frequency use) is 



Fig. 3-10. Normalized loop bandwidth as a function of 
normalized delay 



Fig. 3-11. Technique for measuring natural frequency of a tracking loop 
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frequency-modulated by a sinusoidal modulation m ( t ) = Letting 
A to cos <a m t to produce the signal 

s(t) = A sin [<o c t + 0i (t)] (3-71) 


f * 

Wm “ w w r 2 

r = t 8 /t 2 


where 


and assuming that rr, > > r 2 , we get 


101 (t) = 0 O + f m (z) dz - 00 + — sm <o m (t) (3-72) 

Jo “m 

When this signal is applied to the DAL, the phase detector 
output (the phase error in the Imear model) will be of the 
form 


arg [T (/<4)] = y + tan- 1 <4 


+ . 


x ( , [~ 1 ~ (%/r)sm4'T' T| 
tan - ( w ' mV r) C0S ^rJ/ 


(3-79) 


or equivalently, % is the solution of the transcendental 

Ai.i 

e ( t ) - B m — sin (» m t + 0 m ) (3-73) equation 


The natural frequency is then that frequency < 0 * at which 
the phase detector output is in phase with the input modu- 
lation m ( t ), i.e., d m = tt/2 + 2mr. In mathematical terms, 
we begin by finding the transfer function which relates 
the phase error to the input phase modulation. For the 
linear loop model, 


*(«) = 


A* ( s ) 

1 + AK exp ( — sT a ) F ($)/$ 


(3-74) 


where 


, , On A(l> 


(3-75) 


tan* 




+ u>' m T' — tan- 


J , f 1 ~ (o>m/r) sin »' m T' 

r m ll-(«yr)co S «' m T\ 


| = 2 Ti-n- 
(3-80) 


For T = 0 (zero delay) and t,/t 2 > > 1, i e , 

KM 


tan 1 

the solution to (3-80) is 


0>m 


' -\TF 


(3-81) 


or 

Substituting (3-75) for 0 X ( s ) and (3-46) for F(s) into (3-74) 
gives, upon simplification, 


\T7 


T 2 


*(s) = 7XKTW ( 3 - 76 ) 

S 1 o> m 

where 
T(s) = 

(1/<o ot Ti) [00 (s 2 + dm) + ^A<o] [1 + TjS] 

s 2 + s { [1 + AKt 2 exp ( —sT 8 )]/ n} + (AK/tJ) exp ( ~sT „ ) 

(3-77) 


which checks with a well known result for second-order 
loops (Ref. 3-2, Eq. d -24). For T' > 0 and t,/t 2 > > 1, «4 
will be greater than V"r) and we may still assume that 
tan -1 [<u m (r a /r 2 )] = 7r /2. Using this assumption in (3-80), 
we find, after considerable simplification, that J m must 
satisfy 

1 sin o) m T <o m cos <a m T *1 it 

Lcos J m T + 4 sin <o' m T' - (o/ m 2 /r) J “ T 

(3-82) 


Thus, from (3-76) and the above definition of natural fre- 
quency, we wish to find the frequency <«*, at which 


arg [T (/<,)] = -|- + 2n77 (3-78) 


or equivalently, 

Sin a>mT' — < 0 m cos (%T' < 0 
cos <a' m T' + <o' m sin <o' m T' — (oj^/r) = 0 


(3-83) 
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For small T' (small values of delay relative to r 2 ), the 
solution to (3-83) is 




l 

V (i/r) - v + (ry2) 


(3-84) 


More generally, the solution to (3-83) can be expressed 
implicitly as 



(3-85) 


where x is the inverse solution of 


, x r 

/ 4 / cos x \ 

1 

r = 

. / 1 H ( — ) - 1 

tan x 

2 L 

\ r \sur x ) 

J 


(3‘86) 


Plotted in Fig 3-12 is the ratio of <*' m as given by (3-85) 
to its zero delay value (Eq. 3-81) vs normalized delay 
T' = TJ t 2 , with r as a parameter. As was true in the case 
of loop bandwidth, the presence of delay in the loop can 
significantly increase the loop’s natural frequency 



Fig. 3-12. Normalized loop natural frequency as a 
function of normalized delay 


E. Data-Aided Loop Acquisition 

The problem of characterizing the acquisition proper- 
ties of the DAL bears directly on the loop’s nonlinear sta- 
bility performance discussed previously. In other words, 


in the region of stability, the loop will either acquire 
(achieve a stable lock point) or go into a stable limit cycle 
Unfortunately, because of the extreme difficulties of the 
analysis, we cannot quantitatively characterize the DAL’s 
acquisition performance, however, as before, we can infer 
certain qualitative measures from the acquisition theory 
for PLLs with delay developed in Ref. 3-2. In particular, 
the function 



coso) fT 8 i sina)/T.,"] 

y {0,,) =AK + 2r. 

<j>i T„ oif T ^ J 


(3-87) 


defines the synchronization boundary , where «>/ is the 
fundamental beat note in a Fourier series expansion of 
<t> (t). For any given value of normalized delay 


T 


a T. 

Ti 


setting 

dy M _ q 

3(0 f 


(3-88) 


gives in general several values of <o f , say {w/ n } (see 
Fig. 3-13), which satisfy Eq. (3-88). Calling the smallest 
of these values u>, 0 , and denoting the initial frequency 
offset by Q 0 , the following conclusions may be drawn. 


(1) If fio/AK < y («v 0 ), the loop will lock. 

(2) If Q 0 /AK > y («v 0 ), frequency pushing occurs, and 
the loop goes into a stable limit cycle 


Expanding upon the second conclusion, we note that for 
any given small value of normalized delay T' and a value 
of Oo/AK > u/ o , there will be two solutions to the equation 

~ 7 (<*>/) (3-89) 

namely, y («>/,) and y (w/J, which represent equilibrium 
states with corresponding beat frequencies a >/ t and o>/ 2 . 
The significance of these equilibrium states is that they 
represent the conditions of false lock of the loop. If we 
order the beat frequencies such that w fl > o> f „, then the 
state (d/j which occurs at a negative slope of the synchro- 
nization boundary is unstable, and the state <u/ 2 , which is 
located on a positive slope of the synchronization bound- 
ary, is stable. As the delay in the loop is increased, we see 
that for the same value of loop detuning, more than two 
equilibrium states can exist. These states alternate be- 
tween stable and unstable, much like the lock points of a 
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Fig. 3-13. Synchronization boundary of a second-order 
phase-locked loop with arbitrary delay 


Returning now to the first conclusion, for a PLL with 
small delay T s < < 1, the loop will lock if 


1 tip 

AK 


< 



2rW L T a 
r + 1 


(3-92) 


and F a < < 1. Clearly, when no delay is present in the 
loop, (3-87) reduces to the well known result for the acqui- 
sition range of a PLL, namely, 

-J^L<V2F; (3-93) 

Thus, we observe from (3-92) and (3-93) that the presence 
of delay in the loop reduces the acquisition range 


Experiments have been conducted on the second-order 
sampled data hardware version (Section IV) of the DAL 
to determine the severity and number of false-lock points 
that can occur under a specified set of conditions To this 
end, the extreme case of S = 10 was chosen, together with 
the conditions that r = 2 and AK — I0 4 . Acquisition tests 
were then conducted by offsetting the loop VCO from 
the nominal transmitted frequency by a certain amount 
when the loop was open*, and then closing the loop and 
observing the transition to the steady-state condition 
These tests were conducted in the absence of additive 
noise 


PLL with no delay As the delay is increased still further, 
the loop ultimately will not lock even if the VCO is at its 
quiescent frequency and the offset in frequency is zero 
This corresponds to the situation where the equation 

y (<■>/) — 0 (3-90) 

has at least two solutions. For example, in Fig. 3-13, a 
value of T' = 0.2 represents such a situation. In this 
instance, the loop is said to be absolutely unstable; i.e., it 
can never lock. The critical value of delay at which this 
occurs is that value which simultaneously satisfies the 
equations 

y (o>/) = 0 

(3-91) 

h M _ n 

dwf 


Figure 3-14a shows the results for the normal DAL 
S-curve. The offset frequency Q 0 (rad/s) has been normal- 
ized to the loop bandwidth False-lock points are indi- 
cated by the sohd vertical lines topped by numbers. The 
vertical dashed lines represent instability region bound- 
aries. Some instability boundaries are coincident with the 
false-lock points The arrows indicate the direction taken 
by the loop VCO once the loop is closed, assuming that 
the initial open-loop frequency offset falls somewhere 
between the relevant instability boundaries To illustrate 
this, for all initial | Q„/2ttW l | <0.19, the proper lock point 
is reached. For 0 19< \Q 0 /2irW L | <0 42, the first false- 
lock point becomes the steady-state condition. For 
0.5< J Q 0 /2ttW l | <0.68, the loop will tend to false-lock at 
point 3 The alternate regions of false lock are indicated 
by shaded arrows. The reasons why the false-lock points 
are distributed as shown, or why some instability region 
boundaries are coincident with lock points, is unknown 
Notice the definite subregions where frequency pushing 


This critical value of delay was given in Eq. (3-59) in con- *The condition of open loop involves removing any residual mem- 
nection with our discussion on nonlinear stability. ory from the loop filter F(p) 
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^ — PROPER LOCK POINT 

Fig. 3-14. DAL false-lock ranges for 5 = 10: a. normal S-curve, b. inverted S-curve 


away from the proper lock point occurs. This would be Since the noise-free S-curve of the DAL behaves like 
indicated by any arrows that point to the right sin 2<f> in the vicinity of the origin, and since the S-curve 

ef (<j>) = sgn (cos <f>) appears to give improved acquisition 
A means for minimizing the number of false-lock points range, it is natural to ask if we can make the DAL behave 

for a given 8 and eliminating the frequency pushing is like sgn (cos <j>) in the vicinity of the origin with a simple 

advanced by the following arguments configuration change. Indeed, this is so and is essentially 

accomplished by inserting a minus one gain into the loop. 
The acquisition range for a tracking loop with a rectan- "^ 1 ' s inverts the phase detector characteristic and causes 
gular phase detector characteristic, e.g., of (<j>) = sgn (sin <f>) shift by v/2 rad to achieve stable lock. Thus, the 

and no delay in the loop, requires that resulting noise-free S-curve becomes 

J^L <2 V"f 7 (3-94) ~ sin (* + i) Sgn [ cos ( 0 + T)_ 

= cos (</>) sgn (sin <j>) 

when < < 1 Thus, we conclude from (3-93) and (3-94) 

that a loop with a rectangular characteristic has a wider After the loop has obtained a locked condition [as deter- 

acquisition range by a factor of V" 2 than a PLL with a mined by a lock detector (Section V)], the sign inversion 

sinusoidal characteristic would be removed 


An expression analogous to (3-92) for a loop with a rec- 
tangular S-curve and arbitrary delay Ta is presently not 
available However, from the zero delay results of (3-93) 
and (3-94), we might conjecture that for small delays, 
(3-92) applies if it is multiplied by Y~2. Finally, if the loop 
S-curve is of the form of (<f>) = sin 2 <f> or ©f (<f>) = sgn (cos <f>), 
then the results obtained for the acquisition range of the 
PLL and loop with of (<f>) = sgn (sin <f>) apply, respectively, 
if Cl 0 is replaced by 2f2 0 , i e , the acquisition range of a 
tracking loop whose phase detector characteristic is com- 
pressed by a factor of 2 is reduced by a factor of 1/2. 


Experimental evidence tends to indicate that the in- 
verted and phase-shifted S-curve indeed results in an 
improved acquisition performance 

For the conditions stated above, experiments were con- 
ducted using an inverted S-curve In addition, it was 
found that the open-loop gam had to be reduced by a 
factor of 150, and r 2 increased by a factor of 3 for satis- 
factory performance The results are shown in Fig. 3-14b. 
Note first that the range for which the proper lock point 
is achieved is increased by a factor of about 1 3 Note 
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also that the trend for large offsets is to go in the direction 
toward lock, even though the loop comes to rest at false- 
lock points It has been observed that these false-lock 
points are not very strong, and it is quite easy to make 
the loop leave a false-lock point and seek one to the left. 
This suggests that if a swept acquisition ramp into the 
VCO is used, the loop could be acquired easily. Our final 
note is with respect to the two regions where small, 
shaded arrows point in both directions This is intended 
to convey the fact that in these regions, the loop tends 
to wander between the instability boundary and the false- 
lock point Thus, there appears to be little loop gain, and 
therefore control, in these regions. 

Clearly, very little is known about acquisition for small 
S Determining the best way to achieve lock for large 
frequency offsets in these cases will require further re- 
search and experimentation. 


F. Performance of the DAL When Preceded by 
a Bandpass Limiter 

In the practical operation of a coherent receiver employ- 
ing a carrier-tracking loop, the loop is often preceded by 
a bandpass limiter (BPL) because the loop bandwidth is a 
function of the signal power, and any fluctuation in this 
quantity (owing to a change in range, for instance) causes 
the loop bandwidth to fluctuate similarly. BPLs are also 
used to protect various loop components, the multiplier in 
particular, where signal and noise levels can vary over 
several orders of magnitude and exceed the dynamic range 
of these components In this section, the theoretical results 
needed to explain the behavior of the DAL when pre- 
ceded by a BPL are discussed The detail theory and 
operation of the BPL are covered more completely in 
Refs 3-2 and 3-3. 

A block diagram of the DAL preceded by a BPL is 
shown in Fig 3-15. Since the success of the DAL depends 
on the signal sidebands, the IF filter must be sufficiently 
wide to pass the total modulation. Letting W, denote the 
two-sided bandwidth of this filter, the signal-to-noise ratio 
at the BPL input is 


2 P D 

P ' N 0 W, 


(3-95) 



Fig. 3-15. The DAL preceded by a BPL 


design point loop bandwidth. In the case of suppressed 
carrier transmission, the appropriate expression is 

_ 2 P co 

0 n 0 w L0 

where the zero subscript denotes the value of the param- 
eter at the design point. Defining the BPL suppression 
factor aj by 




(3-96) 


then, using (3-25), the DAL signal-to-noise ratio can be 
expressed as 


P 


= x 



1 + *~o 
1 + r„//x 


) 


where 


A 2P i) erf 2 VlL 
x “ N 0 W lo 


a W L0 
y W, 


(3-97) 


To — «io VP^K erf VRTo^F, 


In designing phase-coherent receivers preceded by 
BPLs, it is convenient to characterize loop performance 
at a given operating point relative to its value at a fixed 
design point. In the past, it has been customary to specify 
this design point by the carrier signal-to-noise ratio in the 


A <*io 



The symbol «i 0 represents the value of <*, at the design 
point, where, by definition, x — Y 0 . Also in the above, 
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Pj = 8/7 r 2 is the total power in the first spectral zone and r„ is defined in Ref. 3-4 by : 


1 ~ exp (~ P% ) 

lp (tt/ 2) (p,/2) exp (- P< ) [fo (p,/2) + 7i ( P ,/2)p {1 + [(4/T*r) - 1] exp [— p« (1 - »/4)]> 


(3-98) 


The p.d.f. of the phase error is still given by (3-16), together with the potential function defined in (3-47). However, the 
parameters a and /? are now expressed in terms of the design point quantities by 


a — 


P = 


fo/p. + 1 \ _ 1 — Fj 

f o/P / P (fo/M) 

r 0 /n + 1 \|~ p / + 1 \~1 

r„/p /(_2W L n\ r 0 /p + 1/J 


[Op- VPoKerfV^Ml 


— Fi) of ($)] + (</>) 


(3-99) 


where p is given by (3-97) Also note that since the BPL appears only in the lower arm of the DAL, the data signal-to- 
noise ratio is unaffected by its presence The entire penalty in data detection performance due to the presence of the 
BPL comes from the increased noisy reference loss. 


*rp is the value of r p at p, = 0, which, for an ideal bandpass filter, is approximately equal to 1/0 862 
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IV. DAR Implementation 


The subsequent sections deal with the implementation 
and performance of the symbol detector (matched filter 
and decision device) and sychronizer, the delay function, 
and the general organization of the data-aided loop 
(DAL) We are concerned with effective, economic imple- 
mentations that will result in a DAR performance reason- 
ably close to that obtained under optimum conditions 

A. Digital vs Analog Mechanizations of the 
Entire DAR 

Before proceeding with implementation details, it is 
necessary to develop a rationale which permits the DAR 
to operate over a large range of conditions such as input 
frequency, data rate, loop bandwidth, etc , without undue 
complications, the need for constant redesign, and with a 
minimum of actual component or block changes Cost and 
reliability are uppermost considerations. 

One approach would make use of modern, operational- 
based analog circuits to achieve the desired functions 
Another would employ as many digital or saturated 
devices as possible, and produce the desired results by 
Computationally operating on samples of the input sig- 
nals)— in effect, a special-purpose computer. The exclu- 
sive use of analog circuits requires a large amount of 
analog signal and component switching when it becomes 
necessary to change conditions. Analog circuits are also 
subject to drift and tolerance problems. On the other 
hand, the digital embodiment of second- and third-order 
loop filters, voltage-controlled oscillators (VCOs), and RF 
phase detectors is quite complex. As a result, the present 
state-of-the-art, coupled with practical experience, indi- 
cates that a hybrid solution is best. 

The phase detectors, which must generally operate at 
RF input frequencies in the range from 1 MHz to several 
tens of megahertz (or even higher for very wideband sys- 
tems), are most easily built in the form of a ring-bridge 
diode configuration. At lower frequencies, it is possible to 
use MOSFET switching elements. The VCO is also con- 
veniently constructed using analog circuitry— usually some 
form of voltage-controlled crystal oscillator. The loop 
filter, whether it be for a second- or third-order loop, con- 
sists of resistors and capacitors, together with some linear 
amplification, and is quite easily mechanized using opera- 
tional configurations. All other functions of the basic DAR 
are readily implemented using logical elements. Fig- 



L I 

Fig. 4-1. Functional DAR implementation 


ure 4-1 shows the DAR functional blocks for a hybrid 
mechanization; digital portions are indicated by broken 
lines. The analog-to-digital and digital-to-analog con- 
verters are labeled, respectively, ADC and DAC. 

B. Matched Filter and Symbol Synchronizer 

When the DAL is locked, the stochastic input to the 
matched filter, (t), consists of the low-pass data symbols 
plus noise. Mathematically, the output of the matched 
filter and decision device configurations for any symbol d t 
may be written as 

I f 1 

d. = sgn < / (t) dt (4-1) 

U 0-1) r I +'« v 

where £ is the symbol synchronizer’s estimate of the true 
symbol epoch e. For convenience, and without any loss of 
generality, we shall assume that e = 0, and thus, e repre- 
sents the relative timing offset due to the symbol synchro- 
nization process. The symbol synchronizer is designed to 
maintain E (e) = 0. 

The digital equivalent of the matched filter and deci- 
sion operations is shown in Fig. 4-2. The low-pass data 
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Fig. 4-2. Digital equivalent of matched filter and decision device 


symbols plus noise are low-pass filtered (specification of 
bandwidth is discussed in Section D) and sampled at the 
Nyquist rate to obtain the samples €1 (t„) The ADC com- 
prises a number of uniform levels L — 2 K , with K the num- 
ber of bits, such that a nearly linear sample estimate of 
the input is made. (See Section E for a discussion of the 
number of ADC levels ) The summer accumulates the 
sample values taken in one symbol period as directed by 
the timing information derived from the symbol synchro- 
nizer. The epoch of the sum is determined by an integer 
/ (e), with E {] (£)} =0 The data estimate is equal to the 
sign of the accumulator at the end of the sum process. 
The accumulator is reset to zero after each symbol. 

The performance of the digital matched filter can be 
made to operate within 0.1 to 0.2 dB of that for an ideal 
matched filter. The major contributors to a few tenths of 
a decibel in performance loss are the jitter on the symbol 
synchronizer timing estimate [i e , the variance of / (£)], 
and any signal-to-noise ratio (SNR) degradation due to the 
finite number of ADC levels Obviously, the closer we 
wish to come to the ideal case, the greater the complexity 
and cost of the mechanization. The problem of symbol 
synchronization is treated first 

The method whereby the phase of the symbol synchro- 
nizer s local oscillator is discretely updated to match that 
of the incoming bit sequence makes use of the digital 
phase-locked loop principle. In particular, the phase detec- 
tor topology of the data-transition tracking loop (Ref. 4-1) 
was found by experimental simulation to provide the best 
error signal, which is used either to bump the phase of the 
VCO clock ±aT from its present position, or to leave 
the VCO at its phase position 

The time increment AT = T,/2N is the time between 
ADC samples. The probabilities with which these three 
possible timing updates occur depend on the magnitude 
of the timing error between the incoming symbol sequence 
and the clock, and the presence or absence of a transition 
in the data. 



Fig. 4-3. Model for data-transition type of digital 
symbol synchronizer 


Referring to Fig. 4-3, the input signal s ( t , e) is a random 
pulse train, representing the data symbols; i.e., 

s ( t , e) = 2 d t u [t — iT s — e] (4-2) 

where u(t ) is the unit rectangular pulse defined by 
u (t) = 1 for 0 ^ t ^ T„, u (t) = 0 for all other t, T, is the 
symbol period; the symbol amplitude d l takes on the 
values ±"V Pd with equal probability, and e is the random 
epoch to be estimated, which, as previously stated, will 
be assumed to equal zero. The input additive noise process 
n(t ) is assumed to be white gaussian, with a two-sided 
spectral density N 0 / 2 W/Hz The in-phase and mid-phase 
filters, F, P (s) and F„ P (s), respectively, are approximated 
by the accumulator shown in Fig. 4-2 and a second 
accumulator which begins its sum of samples in the mid- 
symbol position [i.e , when n = N + 1 + J («)]. The deci- 
sion devices simply extract the sign of the indicated 
accumulator 

The transition detector then examines two adjacent 
decisions, d t . u o„ and records an output 7, according to 
the following rule: 

If d , = d,. u then 7, = 0 

If d, = — 1 ,d t -i = +1, then 7, = +1 (4-3) 

If d x = 4- 1, d ,- 1 = — 1, then 7, = — 1 

The output 7, of the mid-phase filter is also sampled 
at inteiVals of T, and must be delayed by an amount 
r = T,/2 before multiplication with the appropriate 7,. 
The error signal, e t = 7, X J„ is in general accumulated 
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for an interval MT„ (i.e., M values of e t are summed), and 
the clock is bumped (i.e., the phase estimate £ + i at interval 
MT S is updated) according to the following rule: 

e +a =£-AT, if2e, <0 I 

= e if 2 e > = 0 > (4-4) 

/ e +1 = 'e + AT, if 2 e, > 0 ) 

Now e" +1 is a discrete random variable which ranges over 
the set of 2N possible values in the interval T a . If we 
normalize £ +1 by T„ then the normalized random variable 



ranges over the set of integers —N,—N + 1, • , — I, 

1, • , N — - 1, N (see Fig. 4-2) For simplicity of nota- 

tion, we shall subsequently drop the +1 subscript on £,,, 
since it is only relative 

In order to find the performance in the presence of 
noise, it is necessary to determine the steady-state proba- 
bility density function p (J) of J (e) and then average the 
dependence of the symbol probability of error Pj on € over 
this density. The method used is based on a random walk 
model, and details may be found in Ref 4-2. 

The computation of p (}) is rather tedious An example 
is given in Fig 4-4 for M = 1, R s = 3, N = 8 We note 
that E (/} = 0, but that J — 0 is not allowed. As R„-» oo, 



Fig. 4-4. Probability frequency function of J («); M = 1 
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the loop in effect becomes a bang-bang servo, switching 
between / = — 1 and J = + 1 as a result of the frequency 
drift between the transmitter and receiver timing oscilla- 
tors and inherent circuit noise. 

Figure 4-5 shows p (J) for M = 64, with other condi- 
tions remaining the same. Comparing this result with that 
of Fig. 4-4 shows a definite decrease in the variance of /, 
thus indicating that increasing M acts to decrease the 
effective bandwidth of the symbol synchronizer loop. A 
plot of the variance of /, o-J, as a function of R s for the 
cases of M = 1 and M = 64 appears in Fig. 4-6. 

The choice of the size of M in practice depends on 
(1) the desired acquisition and transient response time 
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. 4-5. Probability frequency function of J (e); M = 64 
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(usually measured in symbols) and (2) the allowable deg- 
radation in P® performance relative to ideal, due to a]. 
Tests have shown that 2 < M < 16 appears to give the 
best results from the acquisition standpoint, with M = 4 
being used for most of the system testing described in this 
report. Actually, the choice of M for the acquisition case 
primarily affects the mean time to acquire symbol sync, 
which is generally very much less than the time required 
to obtain carrier synchronization. For M = 4, N = 8, and 
R s = 4, the mean time to acquire symbol synchronization 
has been observed to be about 26 T a . 

Table 4-1 shows a list of measured losses resulting from 
<rj, taken under the stated conditions and inferred from 
the measured P®. Figure 4-7 presents the error probability 
performance as a function of signal-to-noise ratio R s , in 
dB, for the extreme cases of M = 1 and M = 64, obtained 
by averaging P® (/) over the probability density function 
p (/). The error probability performance of an ideal phase- 
shift-keyed (PSK) system is indicated in this figure for 
comparison. Also included is a set of experimental points 
taken for M = 4. Note that they fall closer to the M = 64 
curve than the M = 1 curve, which indicates that rapid 
improvement with respect to M = 1 is possible for small 
M, and that, practically, M on the order of 4 is sufficient 
This conclusion has already been demonstrated by the 
experimental results shown in Table 4-1. 


Table 4-1. Effective matched filter loss as a function of M a 


M 

Measured loss from ideal performance, dB 

1 

1 1 

2 

07 

4 

05 

8 

04 

16 

03 

32 

0 25 

64 

0 25 

*N = 8, R, 

= 4, oscillator offset = 0 01% 


One final point should be made before concluding this 
section. There will normally exist some frequency differ- 
ence between the transmitter and receiver master symbol 
timing oscillators, this has been referred to in Table 4-1 
and Fig. 4-7 as oscillator offset, specified in percent of the 
symbol frequency An offset on the order of 0 0\% has been 
found by measurement to cause a measured P® degrada- 
tion of between 0 1 and 0 2 dB. Some experiments were 
run to determine how large an offset could be tolerated 
while still maintaining DAL lock. For the case of M = 4, 
N = 8, and R s = 4, the upper limit was found to be in the 
neighborhood of 0 5%. 



Fig. 4-7. Theoretical and measured symbol 
synchronizer performance 


C. The Delay Function 

The delay function, which delays the quadrature phase 
detectors output by T„, is conceptually and mathemati- 
cally a very simple device. Its hardware realization, how- 
ever, is another matter. Considering for the moment how it 
might be mechanized without the use of digital com- 
ponents, at very high data rates it could perhaps be 
realized by using a suitable length of wideband coaxial 
line. As the data rate becomes progressively lower, the 
length of the coaxial line will become prohibitively great, 
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req uir ing that some form of lumped element delay line 
equivalent be substituted. Such mechanizations, being 
composed of R, L, and C components, become physically 
large and difficult to construct as the delays approach the 
tens of millisecond region. Another problem is that of 
adjusting the delay lines to give accurate delays, and 
maintaining them over large temperature variations. 
Other types of delay lines, such as magnetostrictive and 
surface wave, could also be used, however, these types are 
not very amenable to large delays (ie., greater than 
100 jus). Thus, the analog delay line is not a very practical 
solution to the delay function problem 

Digital implementations of the delay function are 
straightforward and simple to implement. Consider the 
diagram shown in Fig 4-8. In a manner identical to that 
of the matched filter, the low-pass process is sampled by 
a K-bit ADC, where the sampling clock is identical to that 
used for the matched filter. The same clock is also used to 
operate a 2iV-stage shift register, thus providing the requi- 
site delay. Large or small delays present no problem, since 
they are dependent only on the sampling clock frequency. 

An alternate mechanization is shown in Fig 4-9. The 
2N-stage shift register is replaced by an accumulator 
identical to that used in the matched filter, and with the 
same timing. The output e 3 (t„ — T„) now occurs at a rate 
equal to the symbol rate 0},. Thus, this implementation 
acts to provide the requisite delay, and also to low-pass 
filter the loop phase error-related term The filtering proc- 



SAMPLING 

CLOCK 

2N.» 


Fig. 4-8. Digital equivalent of delay line 
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Fig. 4-9. Alternate digital delay line 


ess has no real effect on the tracking loop performance, 
however, as W L <<£R S for proper loop operation (Sec- 
tion III-D), meaning that the loop filter F (p) determines 
the loop SNR. Measurements have been unable to discern 
any performance difference between this mechanization 
and a true delay line A definite advanatge of the accumu- 
lator mechanization over the shift register is that the 
output signal rate is less by a factor 1/2N, so that the 
subsequent circuit blocks, especially the DAC, need work 
only at the symbol rate We will also see in Section V-B 
that this is convenient to the implementation of the 
lock detector 

D. Choice of Sampling Bandwidth 

Having established that the matched filter and delay 
function should be digitally implemented, we next address 
the problem of sampling bandwidth. As indicated in 
Figs. 4-2 and 4-8, the sampling bandwidth was arbitrarily 
chosen to be N£R S Hence, the question of importance is. 
How small can we choose N without paying an appre- 
ciable penalty in SNR performance'* It is desirable to 
make N as small as practical for the highest symbol rate 
used, as this has a direct effect on minimizing the speed at 
which the ADC has to function, and on the accumulator 
sizes, which in turn minimizes implementation cost. 

In terms of system performance parameters, we aim to 
compute the degradation of SNR with decreasing N and 
determine from this functional relationship whether or not 
a distinct lower bound exists. 

The performance of the matched filter configuration 
shown in Fig 4-2 is the subject of the analysis, where it is 
assumed that the symbol synchronizer acts to provide a 
perfect timing reference. In terms of the input symbol 
stream, we reason that with long strings of ones (or minus 
ones) in the input data, the probability of error of d will 
not be degraded by the finite bandwidth effect of the low- 
pass filter (LPF) Clearly, then, the worst case insofar as 
degradation is concerned will occur for alternating data 
(i.e , 1, — 1, 1, — 1, 1, — 1, ). We shall assume this worst 

case in the analysis which follows by modeling d ( t ) as a 
square wave of period 2T S . 

Consider the system model of Fig 4-10 and its mathe- 
matical equivalent. Fig. 4-11. The signal s\ (#) is a periodic 
impulse train with Fourier transform 

= (4-5) 

k=- <x> 
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This impulse train is passed through a filter whose impulse 
response is 



d(r) 


t 


Fig. 4-10. Sampled data system model 


Taking the Fourier transform of h 0 (t), we have 

H< <fl = ^ sta (l0[ 1 - cos (iO] (4 ' 7) 

The signal component of the output of H 0 (/) is a periodic 
square wave, which meets the requirements on d ( t ). 

Passing s , ( t ) through the filter H 2 (f) provides an im- 
pulse sampling waveform of period 1/(2 N^,). The im- 
pulse response h 2 ( t ) is characterized by 

2 X 

m*) = s ( t ~ 4 n&.) (4_8) 

n = -2N + l 

with Fourier transform 

2N 

<«> 

n = l 

The filter whose bandwidth is to be determined is 

Hi (f) = \ H l (/)| exp [/©(/)] (4-10) 




Fig. 4-11. Sampled data system mathematical equivalent 
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and the output of this filter is denoted by x f ( t ) = s/ ( t ) + 
The sampled filtered square wave appears in 
Fig. 4-11 as Sj (t) and has a Fourier spectrum 

Sx (f) = [S, (f) Ho (f) Hi (f)] * [S, (f) H 2 (f)] (4-11) 


The noise power in the data estimate d can be ex- 
pressed in terms of the variance a\ of the unsampled zero- 
mean noise process at the output of H 1 (f) and the noise 
correlation between samples of this process. 

Defining the normalized filter 


segment of this signal, denoted by s 2 (t), can 

be obtained 

Hi (/) = Hx (N<(R,f) 

(4-17) 

by multiplying s x (() by 



gives 


1 

(t) = < 

1 T 

I 1 ' 1*1 

) T 

(4-12) 

k - 1 


| 

o, |t|> T 


N 0 



( 1 1 2 


0 ti 

~ 2 T s *». 

(4-18) 

Hence, in the frequency domain 


where 


S 2 (f) = 0{s 2 (t)) = S 1 (f)*S 2 (f) 

(4-13) 

ps = 2 | Hi (/) | 2 cos irfk df 

(4-19) 


Since x 2 ( t ) = s 2 ( t ) + n 2 (t) is an impulse sampled version 
of the filtered input square wave plus noise over a single 
symbol interval, integrating x 2 (t) between — oo and oo 
(or, equivalently, between —T ,/ 2 and T,/ 2) is equivalent 
to summing ZN samples of x, (t). Hence, the data estimate 
d may be expressed as 




The signal-to-noise ratio of d is now defined as 



1 / 2n— 1\ 


(4-20) 


where P, and al g are found from (4-16) and (4-18), respec- 
tively Making these substitutions, we find that 


n--N + 1 


= j ~r 

2Nj_, 


x 2 ( t ) dt 


(4-14) 


CMR - /2SWPH 

v No )(.<_ 


(4-21) 


The next step in the procedure is to obtain the signal-to^ 
noise ratio of a. The power in the signal component of d 
can be obtained from Eq. (4-14) as 


The factor P',/ a' n 2 s represents the SNR degradation with N 

Two sample cases have been calculated and plotted to 
show the effects of N on performance. The filter forms 
taken are a cascade of K identical poles, for which 


P. = [E{^] 2 = [^J (4-15) 


|H ' l(f)|2 [i + if/feTY 


(4-22) 


After some straightforward algebraic manipulation, 




2 [(2k-l)/4 N] 


oo 

I 2 = 0 

= aTsYpT 


k + [(2k - 1 )/4N] 


||h^2N£? 8 ( 


u + 


2k — 1 

4N 


[(2k - 1)/4N] 2 } 

= 1 

)]|co,{e[21»«.(l, + ^i)]} 


(4-16) 
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Fig. 4-12. SNR loss as a function of N for a cascade 
of K identical poles 


and a K-pole Butterworth ( K even), where 

l« ; WI‘ = ! + (,//.)«•■ f 4 - 23 ) 

Details pertaining to these filter forms, as well as the 
values of p*, may be found in Ref. 4-3 

Plots of the SNR degradation for the two filter forms 
appear in Figs 4-12 and 4-13 Clearly, the Butterworth 
filter is the best choice. Experimental results for N = 8 
and K = 6 Butterworth filters indicate a loss of about 
0.3 dB. Obviously, from the curves, N should be taken to 
be 4 or greater in order to keep the losses at a reasonable 
level. The implementation descnbed in Section F takes 
N = 8. 

E. ADC and DAC Sizes 

The process of analog-to-digital conversion is referred 
to as quantizing of the input signal plus noise ensemble. 
We wish here to focus our attention on the number of bits 
required by the ADCs and the DAC. 



Fig. 4-13. SNR loss as a function of N for a K-pole Butterworth 

Two questions arise relative to the quantizing operation. 

(1) How many bits K should be used m order that the 
SNR of the quantized samples be insignificantly 
degraded relative to ideal samples' 1 This loss will 
be referred to as the quantization loss 

(2) What is the effect of the discrete nature of the 
matched filter accumulator on the probability of 
symbol error, as a function of K ? 

The answer to these questions is important, as it essen- 
tially determines the word size of the accumulator opera- 
tions which follow the ADC Quite obviously, the fewer 
bits per sample, the less the number of total bits of stor- 
age required. However, if one lets K = 0 (quantizing the 
sign of the sample only), it can be shown that for low 
SNRs (and gaussian noise), the SNR is degraded by about 
2 dB, an obviously excessive loss We thus must determine 
a reasonable value for K. 

We denote the output of the LPF preceding the ADC 
as X ( t ) Since the process of quantizing can be viewed as 


JPL TECHNICAL REPORT 32-1583 


35 





a zero memory transformation of the continuous variable 
X ( t ), the processes of sampling and quantizing can be 
interchanged, as shown in Fig. 4-14 (i.e., we let N.£R a = W). 
Further, for the purpose of analysis, since the sampling is 
taken to be at the Nyquist rate so that the samples repre- 
sent the continuous process contained in the band (0, W), 
it is expedient to work with the continuous quantized 
process. Thus, we consider the model shown in Fig. 4-15. 



Fig. 4-14. Equivalent quantizers 



Fig. 4-15. Quantizer analysis model 


Now the quantizer Q is represented by the transforma- 
tion characteristic shown in Fig. 4-16, where the staircase 
is ideally assumed to have an infinite number of steps. 
The finite bit quantizer then becomes approximated by 
the ideal situation when the peak signal plus noise values 
at the input are constrained (scaled) to be less than the 
2 K th level. 


By expanding the derivative of the quantizer character- 
istic in a Fourier series, it can be easily shown that the 
desired input/output relationship for Fig. 4-15 is 



(4-24) 


where e* is the Neumann factor. 

Now, X ( t ) can be written approximately as 

X(t) = VYd ( t ) + n, ( t ) (4-25) 

where the variance of n f ( t ) denoted by var {n f (t ) } equals 
N 0 W. We have neglected the fact that the LPF preceding 



Fig. 4-16. Quantizer characteristic 


the ADC filters out some of the high-frequency compo- 
nents of d ( t ), thus causing a finite transition time between 
symbols. To do otherwise would make the analytical solu- 
tion highly intractable. 

Our goal is to determine the minimum number of 
required ADC bits (K) that produce a quantization loss 
on the order of 0.3 dB Now, if only one ADC level is used 
(i.e., K — 0), such that |X (f) | < V for all t, we find that 

Xk (t) = X 0 (t) = A sgn [X (t)] 

Under this condition, as already mentioned, the SNR loss 
is 2/ 7r or 2 dB (Ref 4-4). Thus, we must have a sufficient 
number of levels to accommodate | X (t) \ max > > V. We 
may then write 

fc = 0 

A V'' €* . r27rkVS' _ . . 2iTrk . . 1 

= 27Z-, I sm [ V + ~V~ n r Wj 

k = o 

(4-26) 

Expanding the sin [ ] term, noting that d ( t ) can take on 
only + 1 or — 1 values, we obtain 

Xk (f ) = TS [ d ^ sin m cos (it n < 

k-0 

+ cos(^^)sin (^- «,(*))] (4-27) 
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Separating the k = 0 term from the expression, we have 
X*(f)=y liTd(t) + n,(t)] 


We must remember that Eq. (4-32) is valid for, say, 
K — 2. Now, a quick examination of the exponential term 
for K = 2 shows that it can be neglected even for k = 1, 
thus, 



E ir . /2 t rkMS\ / 2 A 

J [ d <*) sm y “ J cos n, 



var 


( 7T k U / _ 2ir 2 fc 2 12 


(4-33) 


+ cos 


( 


27rkV^ 

V 



2Jc 

V 



(4-28) 


a well known result for the variance of the quantization 
error 


where we recognize the first term as being nothing more 
than the input function scaled by A/V, leaving the sum- 
mation to represent all of the additional terms which 
degrade the SNR of the samples. We must therefore ob- 
tain the variance of the summation. Evaluating the neces- 
sary integral ( n / assumed to be gaussian), we obtain 


k=i 

x {!-^^H-(^)*T]} <«■> 

Before proceeding with the final step, it is necessary 
to relate the peak value of | X (t) | to the maximum level 
taken for the ADC. For a K-bit ADC, the maximum quan- 
tization level is 2 K V. Now, if p AD c is defined as the input 
SNR to the ADC, then assuming that p ADC < 1, the peak 
value of | X (t) | will be on the order of the peak value of 
the noise, which we will take to be 


4 a n , = 4 i~WoW = 2 K V (4-30) 


Thus, we obtain 


4 V N 0 W 

2 K 


(4-31) 


which, when substituted into Eq (4-29), gives 


Jc = l k = i 

X 1 1 - cos (2*7 rk r^o) exp[-Y (2*-k) 2 ]j (4-32) 


Evaluating the signal mean and noise variance from 
Eq. (4-28), the final result for the sample SNR, defined by 
the symbol p K , is 



pADC A T 

P* — J 2-2JT — pADC I^ADC 

(4-34) 

Table 4-2 gives several values of L AD c- 



Table 4-2. ADC SNR losses 


K 

2 3 4 

5 

L A dc, dB 

-0 26 - 0.06 - 0.02 

-0 01 


It can be seen from the table that if K > 3, the SNR 
loss is negligible. 


We now turn our attention to the discrete nature of the 
accumulator distribution and its effect on symbol prob- 
ability of error. Unfortunately, no simple analytical solu- 
tion for this problem has been readily obtained. Consider 
the following dilemma For a K-bit ADC, the discrete 
probability density function (p.d.f.) of the S + N samples 
consists of L = 2 K+1 — 1 discrete values. If M + 1 sam- 
ples are accumulated, then the resulting discrete p.d.f. of 
the accumulator has 2 U L — ( 2 U — 1) discrete values. The 
actual numerical magnitude of the accumulator p.d.f. 
values depends on the M-fold convolution of the input 
p.d.f , or alternatively, on raising the generating function 
to the Mth power. In either case, the procedure is very 
complex, as it involves a large number of combinations 
for realistic K and M. For example, if K = 5 and M = 15, 
which is typical for an actual mechanization, the number 
of possible accumulator values is slightly in excess of 
2 X 10 6 Actually, it is not necessary to calculate this many 
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values, as it can be shown that only a few thousand are 
really significant to the desired outcome, however, even 
this is an arduous job As a result, the answer to the prob- 
lem has been obtained empirically. 


magnitude plus sign) devices with a throughput rate of 
1.25 MHz, providing for a symbol rate as high as 78 kilo- 
symbols per second (ksps) 


Table 4-3 presents some representative R s losses in- 
ferred from measured P® for specified values of R„ when 
M = 15 We note that the largest losses occur for the 
higher values of fl 8 . This is not at all surprising, since the 
probability of error is determined by the extremes of 
the distribution of the noise, which is significantly modi- 
fied by the fact that beyond some ka n (say, k= 4), all of the 
values of the noise are mapped into the largest quantized 
value. The only way to reduce the loss for small values 
of K is to increase M (and therefore increase the LPF 
bandwidth and the sampling rate). This has been expe- 
rimentally proven for the digital command detector 
(Ref. 4-5), where K — 3 and M = 127, and the loss is on 
the order of — 0 1 dB 


Table 4-3. Measured R, losses inferred from P'l 


\R.,dB 

K 

40 

m 

60 

70 

80 


4 

-06 

' / 1 y ’ | [ 

-06 

-07 

-08 


5 

< -0 1 


- 01 

-03 


6 

< -0 1 


7 

< -0 1 


Clearly, if K==7, while M = 15, no significant degrada- 
tion will result. In the implementation described in Sec- 
tion F, K was taken as 7. 


The function of the DAC is to reconstruct an analog 
form of the loop error signal No studies were actually 
conducted to size the DAC Since K = 7 was the final 
choice for the ADC size, and 2 N — 16, the DAC size was 
simply chosen to be 11 bits plus sign. 


F. Complete Mechanization 

The digital portions of the data-aided loop have been 
implemented according to the block diagram shown in 
Fig. 4-17. As can be seen, two ADCs were used, one in 
each data stream, rather than a single ADC multiplexed 
between the two low-pass phase detector outputs. This 
allows a lower master clock rate for a given symbol rate, 
thus lowering the duty cycle of the ADCs. The ADCs 
selected for use in the data-aided loop are 8-bit (7-bit 


The master clock frequency is 160 times the symbol 
rate This frequency is passed through a series of dividers 
and counters, yielding two pulse trains. The pulse trains 
are present in different phases of a frequency equal to 
1/10 the master clock rate, and hence, cycle 16 times per 
symbol period. 


One of the phases is passed through a 16-bit counter, 
and the occurrence of a 0 or 8 count is detected By logi- 
cally ANDing these detected counts with the generated 
pulse trains, different phases and subfrequencies of the 
master clock are generated. These are necessary to imple- 
ment the digital integrate and dump circuits. Three inte- 
grate and dump, or accumulator, circuits are used. They 
are the mid-phase and in-phase for the matched filter and 
the one-symbol delay. Each of the accumulator circuits 
consists of a 12-bit binary adder and a buffer The accu- 
mulators sum 16 ADC samples per symbol period and are 
cleared by the appropriate clock pulses The mid-phase 
accumulator is cleared by the 8-count detected bit, and 
both the in-phase and 1-bit delay accumulators are cleared 
by the 0-count detected bit. 


The sign bit of the in-phase accumulator is the data 
estimate At each occurrence of the 0-count bit, the new 
data estimate is compared with the previous value. If a 
change is detected, the mid-phase accumulator is multi- 
plied by the data estimate, and the result is added to the 
clock bump accumulator After every fourth change is 
detected, the clock is mcremented or decremented 1/16 
of the symbol period, depending upon the sign of the 
clock bump accumulator, which is subsequently cleared 
Thus, symbol synchronization is maintained. 


The output of the one-symbol delay accumulator is 
multiplied (complemented) by the data estimate, with 
the result being input to a 12-bit (11 bits magnitude plus 
sign) DAC. The output of this DAC is the loop error signal 
which is input to the loop filter. 


The digital logic of the data-aided loop was built en- 
tirely of 14- and 16-pin dip integrated circuits (Series 54, 
74) which were assembled on three (64-socket) cards, as 
shown in Fig. 4-18. The power supplies for the logic and 
converters are external. 
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Fig. 4-17. Complete DAR logic mechanization 
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G. Performance Test Results 

In the preceding sections, we have discussed details of 
the DAR implementation, and presented experimental 
results as they pertained to the immediate subject under 
study. We present here the overall test results as obtained 
from measuring the symbol probability of error The 
experimental data are shown in Fig. 4-19. The tests were 
run at a symbol rate £R S = 1 ksps using a pseudo-random, 
bi-orthogonal data stream The DAL tracking bandwidth 
was chosen to be 20 Hz. 

The experimental data points were obtained for three 
sets of conditions. The first set involves “hardwired” RF 
and symbol sync references, and the corresponding data 
points (indicated by squares) are therefore a measure of 
the LPF and ADC losses As can be seen, the average 
loss attributable to this condition is about —0.3 dB. 

The second set of experimental points (indicated by 
hexagons) was obtained for hardwired RF references, 
but with the symbol synchronizer operating (M — 4) and 
the symbol rate clocks offset by 0 01$. The data points 
average about 0.8 dB offset from the ideal PSK perform- 
ance curve, indicating that the performance degradation 
introduced by the noisy symbol clock is —0.5 dB. This, of 
course, could be decreased by 0 2 to 0.3 dB by increas- 
ing M (Section B). 

The third set of data points (indicated by circles) is 
taken for the complete operating DAR, for which the RF 
references are also noisy. As can be seen, the additional 
loss is between —01 and —02 dB, as the measured rms 
phase error ranged between 4 4 and 8.5 deg. 

The worst-case signal-to-noise ratio m the loop track- 
ing bandwidth occurs for R, = 4, and is 20 8 dB. All of 



Fig. 4-19. DAR performance with bi-orthogonal data 


the data points were obtained for zero detuning, a loop 
gain K = 50,000, and a loop damping factor of V~2/2. No 
limiter was used, as per Fig. 3-14. Additional operating 
conditions are shown in Fig 4-19. 
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V. Attendant Functions 


Two important functions usually associated with co- 
herent receivers are (1) automatic gain control (AGC) 
and (2) in-lock/out-of-lock detection In the following, 
we discuss the implementation and performance of these 
functions with respect to the data-aided receiver 

A. Automatic Gain Control 

In discrete carrier receivers, the AGC voltage is de- 
rived from an in-phase detector, which in effect gives a 
voltage proportional to the carrier level (see Refs 5-1 
and 5-2). The phase reference for this detector is shifted 
90 deg relative to that of the carrier- tracking loop. A sim- 
plified block diagram of an AGC control loop is shown 
in Fig. 5-1, where we have assumed a perfect carrier 
reference. 


GAIN-CONTROLLABLE 



POINT 

BIAS 

Fig. 5-1. Simple AGC loop for discrete carrier receivers 


put of the in-phase detector for the data-aided receiver 
(DAR) is given in Eq (3-8) and is proportional to 
V Pd d ( t ) cos </> ( t ) + N i [ t , </> (#)] If d ( t ) is random (as as- 
sumed), with an equal number of +1 and —1 symbols, 
then the average value of the in-phase detector output is 
zero. Thus, no AGC voltage is generated at this point. 

To circumvent this problem, we deal with the abso- 
lute value of the phase detector output. Because of the 
digital mechanization discussed in Section IV-B, it is 
more convenient to deal with the absolute value of the 
data symbol accumulator (matched filter) output, which 
is obtained by ignoring the sign bit. Thus, the basic AGC 
voltage signal V will be defined as the absolute value of 
the integrated data symbols, whose mean is defined by 

/• V'*' , 

ju.v = E {V} — E I {\Pod (f) cos $ (t) 

+ NAt,<f>(t)]}dt | (5-1) 

For simplicity, we deal here with continuous rather than 
discrete variables and also assume the first transmitted 
symbol. 

For the case of large loop signal-to-noise ratio (SNR), 
we can assume perfect carrier and symbol synchroniza- 
tion, i e , e = 0, cos = 1. Thus, we have 

fxv = E I [ ’ {yfPo d (t) + n 2 (t)} dt \ (5-2) 


The purpose of the AGC loop is to maintain a certain 
degree of linearity through portions of the receiver, espe- 
cially those that handle the data components of the re- 
ceived signal Since the received signal level may vary 
100 dB or more throughout the course of using the re- 
ceiver, the AGC loop prevents the high-gain RF and IF 
amplifiers from saturating (clipping) when the received 
signals are very strong. 

We see from Fig. 5-1 that for large R a , the open-loop 
AGC voltage derived from a discrete carrier is (in the 
absence of a noisy phase reference) a nearly linear func- 
tion of the carrier voltage For the suppressed carrier 
signal, however, there is no carrier voltage. The out- 


where n 2 ( t ) is low-pass white gaussian noise with single- 
sided spectral density, N 0 W/Hz. Now, the integral is a 
gaussian random variable, with mean V Pd T s d ( t ), where 
d (t) is either +1 or —1, and variance N„T S / 2. Refer- 
ence 5-3 shows that 



We see that fi v is a function of both the si gnal and 
noise voltages, and that its minimum value, i e , V N 0 T s /ir, 
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occurs for R s = 0 (no signal, noise only). Thus, it is more 
convenient to examine the normalized mean, viz., 


M'v = YW^/7 = ex P( _R ») + V^erf (\[R.) 

(5-4) 

A plot of Eq. (5-4) is shown in Fig. 5-2. As can be seen, 
the AGC voltage is a nonlinear function of R s . The curve 
becomes quasi-linear for R s > 6, where we can write 




exp ( - R,) 
2R g 


(5-5) 


and for R g > 10, 


r^R, (5-6) 

Thus, under very strong signal conditions, the AGC loop 
for the DAR operates no differently than its discrete 
carrier counterpart. For very weak signals, the AGC volt- 
age becomes limited by the receiver noise. This latter 
condition is actually an advantageous one compared to 
the threshold behavior of the discrete carrier AGC loop, 
where the mean AGC voltage is a strong, nonlinear func- 
tion of the tracking loop phase noise For the DAR AGC 
case, where £ and cos <j> ( t ) become significant factors 
because P„ < < o£ 2 , the dependence of fi v on these factors 
is very small, and thus we still have AGC control even 
though the tracking loop is out of lock. 

It is likely that more filtering of the AGC voltage than 
that provided by the matched filter integrator in the DAR 
is desirable. This can be easily mechanized by adding an 
additional accumulator to the magnitude output of the 
in-phase accumulator buffer of Fig. 4-17 The output of 



R s (NUMBERS IN PARENTHESES ARE IN dB) 

Fig. 5-2. Normalized DAR AGC voltage as a function of R, 


this accumulator may then be converted by a digital-to- 
analog converter (DAC) to provide the gain control volt- 
age. If very large gain variations are desired, the most 
significant bits (MSBs) of the AGC accumulator may be 
employed to switch in or out fixed attenuator pads in 
the RF and IF amplifiers, while the less significant bits 
(LSBs) provide, through a DAC, proportional control in 
the usual manner. Figure 5-3 shows a functional diagram 
of the DAR AGC loop 


MATCHED FILTER 



Fig. 5-3. DAR AGC loop 


B. Lock Detection 

To design a lock detector that indicates when the DAL 
is in or out of a locked condition, one must derive from 
the loop a voltage which on the average is large when the 
loop is in lock and commensurately small when the loop 
is out of lock Then, by comparing this voltage to a preset 
threshold (which may be set based upon the probability 
of detecting the in-lock or out-of-lock state), one may 
detect the presence or absence of lock. 

One convenient method for deriving the lock indica- 
tion voltage is to make use of the AGC voltage V, whose 
mean is defined in Eq (5-1), and its quadrature counter- 
part Vq, which may be derived from the symbol delay 
accumulator (Fig 4-17) and whose mean is defined by 


fj-v Q — E {¥<>} 


jT’ +£ V ypPo d ( t ) sin 4 , (t) + N 2 [t, * (*)] } dt 


(5-7) 


By accumulating V and V Q over a large number of sym- 
bol intervals, and then taking the difference of these two 
quantities, one obtains a voltage which satisfies the above 
requirements on a lock indicator. A lock detector designed 
on this principle is illustrated in Fig. 5-4. 
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SIGN 



Fig. 5-4. Lock detector 


The performance of this configuration is analyzed as 
follows If we assume that when the DAL is out of lock, 
the phase error varies such that it can be assumed to be 
constant over a symbol time, then, from Eqs (5-1), (5-3), 
and (5-7), we have 


the mean of V A is 

m v * = \ { ex P ( - cos2 4>) 

+ exp — R s cos 2 <j> ^ 1 

— exp ( — R s sin 2 4>) 

— exp £ — R a sin 2 <f> 

+ Y V ttR., | cos erf (yfW, cos <f>) 

+ <* JS , t (l-^)erf[vTCco S *(l-^)] 

— sin <j> erf sin <j>) 



Mv = Y^/ _ ~^{ ex P(“ Ra ^ & 

+ exp -R, cos 2 4> (l - 
+ Y VTo T s cos <t> I erf ( V R, cos <£) 

+ ( 1 ~ ^tT) Grf [ ^°° S * ( 1 “ ^)] } 

(5-8) 

MVp = ex P ( “ sin2 4>) 

+ exp[-R 8 sin 2 ,(l-M)]} 

+ ~ yrd sin <t> | er f (Vr; sin 4>) 

+ i 1 - erf [VH7sm*(i - ^f)]} 


- sin * ( 1 - ^)erf [ sin </> (l - ^)]} 

(5-10) 

When the loop is out of lock, we may assume that the 
phase error <j> is uniformly distributed, thus, from (5-10), 

m v * = 0 (5-11) 

v d 

Now, when the loop is locked, we may assume, as was 
done in the AGC case, that for large loop SNR, e = <f> = 0. 
Thus, 

M v * = M a [exp(— R,) + V^erf (Y7Q] (5-12) 

v d 

When the number of accumulation symbol intervals 
M a is large— i.e , greater than 10 (as would be the case in 
lock detection for a carrier regeneration loop)— the nor- 
malized voltage VJ is approximately gaussian distributed 
with the mean given by (5-10) and variance 


Accumulating V and V e over M A symbol intervals, and 
defining the normalized difference variable 


c r w = 


M A 17 

N n T a 


fa 2 

L v 


+ a 2 
v v n 



U A U A 

2V- 

t = l t = i 

V No T f A 


(5-9) 


= -^ {E {V 2 } + E {V* 9 } - [„ v y - [^] 2 } 

(5-13) 


JPL TECHNICAL REPORT 32-1583 


45 








where 


E {V 2 } 
E{V$} 


P d T% cos 2 <f> 

_ 1 + ( 

2|t|' 

V'l 

_L_ 

2 

< Ts , 

) J 

1 

P d T; sin 2 <£ 
2 

_ 1 + < 

(, 2 13’ 

k Ts 

)■] 

+ 


(5-14) 


Again, when the loop is locked and the loop signal-to- 
noise ratio is large, we may assume that <£ = "£ = 0, m 
which case (5-13) reduces to 

<r 2 v N = M a {tt + ttR,- [exp ( — R s ) + V ttRj erf \Tfis] 2 — 1} 

(5-15) 


For very large R„, Eqs (5-12) and (5-15), respectively, 
simplify to 


that the loop is in lock is much much greater than the 
probability that it is out of lock, Eq. (5-18) can be approxi- 
mated by 


P.= prob{v;>V' , i 

V in lock * 



or 


(5-19) 


where n v N and <t 2 n are given by (5-12) and (5-15), 
respectively. 


1 1 r 

— erf 

r^ = v 'n 

d 


9 9 




r 

a. 

1 


and 


= m a V uR s 

(5-16) 

1 

5? 

il 

(5-17) 


If the normalized difference variable V'J is compared 
to a threshold V' and we decide that the loop is in lock 
when V N d > V' and out of lock when V d <V\ then the 
probability of correctly detecting the loop’s state is 
given by 


P. = p,ob{ vj > V' | x ,} X prob{!“ P o “) 


+ prob | V * < V' 


l 


j loop is | 


oToMou. ( X Pr ° b l out of lock / 

(5-18) 


If the threshold is chosen between (perhaps midway) 
the values of fx v n when the loop is in lock and when it is 

out of lock, as given by Eqs (5-12) and (5-11), respec- 
tively, then the two conditional probabilities of (5-18) will 
be of the same order of magnitude Since the probability 


For large R„, (5-19) simplifies to 



- m a - v' ~ 

. V 2M-i (t7 — 1) _ 


(5-20) 


Thus, for given values of R, and P c , and a specification on 
the time to detect the in- or out-of-lock state of the loop 
(equal to M A T S ), one can specify the parameters M A and 
V' For example, consider as suggested that V' — p- v n/2 

and P; = 1X 10 2 . Then, from (5-19), it is calculated that 
for M a = 10, P r = 0 999989 

In conclusion, it should be noted that P r is not a func- 
tion of the absolute signal and noise levels, this is typical 
of a lock detector of the type illustrated m Fig. 5-4. The 
performance is only a function of the choice of V' rela- 
tive to the value of fx >■ corresponding to the worst-case 

v d 

value of R„ and the accumulator size M A Increasing M A 
clearly improves the probability of correctly determining 
the loop’s lock state at the expense of reaction time in the 
determination of either state 
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VI. Complete System Configurations 


In the following sections, we present outline discus- 
sions on how the data-aided loop (DAL) can be adapted 
to serve a number of specific tracking and demodulation 
functions. The intent here is only to define how the loop 
may be configured and to point out some of the salient 
features. The reader may then build upon these basic 
uses to obtain systems that satisfy a particular need. 

A. Modifications of Existing Discrete Carrier 
Tracking Receivers to Achieve DAL 
Operation 

Existing transponder and ground discrete carrier- 
tracking receivers may be easily modified so that they be- 
come suppressed carrier-tracking receivers of the DAL 
type. To illustrate this, a typical spacecraft transponder 
receiver is shown in Fig. 6-1, where the additions are 
drawn as dashed lines and existing blocks that may have 
to be modified are shaded. A ground receiver would be 
modified in a similar fashion. 

The DAL implementation block embodies the functions 
shown in Figs. 4-1, 4-17, and 5-3. One existing receiver 
block that definitely would have to be changed is the IF 
bandpass filter (BPF). In most receivers this filter is quite 
narrow, and would have to be replaced by one whose 
bandwidth is sufficiently wide to pass the data modulation 
sidebands. It may further be necessary to modify or re- 


place the IF amplifier following the filter, as this may also 
have undesirable narrowband characteristics. The auto- 
matic gain control (AGC) filter may also have to be 
changed, so that the AGC loop discussed in Section V-A 
will function properly. Finally, if a “tight” AGC loop is 
designed, it may be possible to eliminate the limiter 
entirely. If the limiter is retained, however, it should also 
be modified to eliminate any narrowband properties. 

For an uplink DAL, an interesting possibility exists of 
resolving the carrier phase ambiguity (mentioned in Sec- 
tions II-C and III-C) during initial acquisition Assuming 
that a lock detector of the type described in Section V-B 
is implemented, it can be used, when the loop is out of 
lock (e.g., during acquisition), to disable the data symbol 
estimate^from the formation of the loop error function by 
forcing a equal to one (Fig. 3-2). Under these conditions, 
the loop is a normal discrete carrier-tracking loop, with 
the exception of the delay element T„ To acquire the 
uplink, then, an unmodulated carrier is transmitted, to 
which the loop locks at one unique point in the carrier 
cycle Now, when the lock detector indicates the m-lock 
condition, the data symbol estimate is enabled, and the 
loop functions as a DAL, with the transmitted data sym- 
bols being a series of ones The loop switches its effec- 
tive S-curve from sin to of (</>), as given by Eq. (3-28), 
without perturbing the lock point. Later, when uplink 
information is to be transmitted, the carrier reference is 
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at the proper phase to demodulate the symbols with- 
out ambiguity. 

B. Subcarrier Tracking 

The use of the DAL is, of course, not constrained to 
carrier frequencies. It may also be used to track sup- 
pressed subcarriers that are encountered in discrete 


carrier systems, and suppressed carrier systems having 
more than one data channel (Section C). Assuming that a 
discrete carrier is modulated by a subcarrier which is, in 
turn, biphase modulated by the data, the configuration 
shown in Fig. 6-2 can be used to track and demodulate a 
sinusoidal subcarrier. Figure 6-3 indicates how a square- 
wave subcarrier might be handled. The only real differ- 
ence between these two configurations is that the type 2 


DISCRETE CARRIER TRACKING LOOP 



Fig. 6-2. Type 1 subcarrier DAL 


DISCRETE CARRIER TRACKING LOOP BPL 



Fig. 6-3. Type 2 subcarrier DAL 
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subcarrier DAL (Fig. 6-3) is easily able to recover the 
harmonic power of the square-wave subcarrier without 
undue specifications on the BPFs, while the type 1 sub- 
carrier DAL, if used with square-wave subcarriers, either 
filters out all of the harmonic power or, if the subcarrier 
BPF is removed, has to operate with extreme limiter sup- 
pression factors 

C. Two-Channel Configurations 

Equation (2-8) gives the general expression for two- 
channel suppressed carrier modulation. The situation for 
quadriphase, with two equal power, equal data-rate 
carrier orthogonal modulations, is a special case and is 
covered in Section D We are concerned here with the 
two-channel cases in which there is a significant differ- 
ence between the individual modulating data rates 

A useful signal form for keeping the two data channels 
spectrally separate is 

s ( t ) = V 2P m di (<) cos <a c t + V 2 P D2 d 2 ( t ) 5 ( t ) sin <o c t 

( 6 - 1 ) 

which may be obtained from (2-8) by letting 

5 1 (t)= 1 

and 

5 2 (t) = 5 (t) 

In (6-1), the data modulation d x ( t ) has symbol rate $R S ,, 
d 2 ( t ) has rate £7? 82 , and we assume that cR ei > > £R a2 , 
implying that P D1 > > P D2 . The square-wave subcarrier 
5 (t), onto which d 2 (t) is biphase modulated, is selected 
so that there is as little mutual interference between the 
two data signals as possible within the overall bandwidth 
constraint placed on the system A typical spectrum for 
s ( t ) appears in Fig. 6-4 



Fig. 6-5. Two-channel suppressed carrier functional 
data'-aided receiver 


With this signal format, the receiver system is essen- 
tially the same as that shown in Fig. 6-3, with the excep- 
tion that the discrete carrier-tracking loop is replaced 
with a DAL which tracks the high-data-rate signal com- 
ponent. A functional diagram is shown in Fig. 6-5. The 
combination of bandpass and lowpass filters in each loop 
configuration effectively eliminates mutual interference. 

There are situations in which it may be advantageous 
always to have the low-power, low-data-rate signal 
directly modulate the carrier in the two-channel mode. 
Consider as an example an uplink to the spacecraft. At 
times, this link is operated single-channel, i.e., command 
modulation only, while on other occasions it is operated 
in a dual-channel mode, i e., command plus ranging 
modulation 

Typically, the command component has considerably 
less power than the ranging component. Thus, the 
superior data-aided receiver (DAR) operating conditions 
would be obtained by tracking the ranging signal. This 
would mean that the ranging signal would directly modu- 
late the earner, while the command data would have to be 
modulated onto an appropriate subcarrier. This, however, 
presents certain receiver reconfiguration problems when 



Fig. 6-4. Typical two-channel suppressed carrier spectrum 
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changing from single-channel (command only) to dual- 
channel (command plus ranging) operation, and vice- 
versa To change from single- to dual-channel operation, 
the receiver would have to recognize that the signals had 
changed, and then switch the carrier regenerating DAL 
to the ranging symbol rate, reacquire, and allow the 
subcarrier regenerating DAL to take up the command 
demodulation/detection function. Not only is this a'com- 
plex switching situation to both recognize and implement 
in the receiver, but it also means that uplink carrier lock 
has to be broken for a period of time, thus disrupting two- 
way tracking 

The alternate, and perhaps more realistic, solution is 
to continue transmitting the command data directly on 
the carrier, and place the ranging on the subcarner. Thus, 
even though the DAL signal-to-noise ratio (SNR) is then 
reduced by the amount of transmitter power allocated to 
the ranging, the loop does not require reconfiguration, and 
does not lose lock. If, in the two-channel case, the DAL 
SNR falls below some desired minimum, it is possible to 
form a hybrid DAL wherein the ranging subcarrier DAL 
may be switched into the command DAL by the subcar- 
rier DAL lock detector to aid and thus increase the SNR 
of the carrier regenerating loop. This concept is discussed 
in Section III-A. 

D, Quadriphase and Higher-Order Tracking 
Loops 

The extension of biphase modulation of a carrier to 
more than two phases is referred to as polyphase or 
N-phase modulation, where N is the total number of 
transmitted phases. In particular, during a given trans- 
mission interval of T seconds, the transmitted signal can 
be expressed in the form 

s(t) = V~2PnSin ^ ^ " j ; 

* = 0,1,2, ,N — 1 (6-2) 

where u> c is again the carrier radian frequency. For N = 2, 
the above signaling format represents phase-shift-keying, 
as discussed in Section II, while for N = 4, it corresponds 
to quadriphase-shift-keying Actually, for any value of 
N > 2, we refer to the signaling technique as multiple- 
phase-shift-keying (MPSK). In the quadriphase case, the 
transmitted signal in (6-2) assumes the form 

s ( t ) = V Pd [d-i ( t } sin a > c t + d; ( t ) cos a> c t ] (6-3) 
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where d,(t) and d 2 (t) are ±1 digital waveforms whose 
transitions occur at T-second intervals For quadriphase 
signaling, the above equation suggests the modulator 
depicted in Fig. 6-6a Setting N — 8 (i.e., octaphase mod- 
ulation), it is easy to show that the circuit in Fig. 6-6b 
generates an octaphase signal. In this figure, d t (t) and 
d, ( t ) also correspond to data sequences of plus or minus 
ones. The generalization of the transmitter modulator to 
N phases is straightforward. 

If one assumes that the channel adds white gaussian 
noise n(t) of single-sided spectral density N„ W/Hz and 
a phase shift to the signal s(t), then the received signal 
can be characterized by 

* (t) = V 2Pz> Sin 4> (t) + - 2 -~ ^ - j + n (t) (6-4) 

where <I> (t) = u>< ; t + 6„, with 8 n a uniformly distributed 
random phase. 

The reconstruction of a carrier from x ( t ) can be accom- 
plished in many ways. For example, the squaring and 



d 4 (0 

Fig. 6-6. Modulator mechanization: a. quadriphase, 
b. octaphase 
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Costas loops used for tracking a biphase-modulated 
carrier can be generalized to Nth power loops and 
N-phase Costas loops, respectively (Ref 6-1). We are 
interested here, however, in the generalization of the 
biphase suppressed carrier DAL to an N-phase configu- 
ration. In what follows, we shall demonstrate one con- 
figuration which is directly applicable to quadriphase 
signaling and then a more general configuration that is 
useful for all values of N 

The biphase suppressed carrier DAL can easily be 
modified to accommodate a quadriphase input signal 
Since, from (6-3), the quadriphase signal input can be 
modeled as the sum of two orthogonal, equal power sup- 
pressed carrier signals, one might anticipate that the nec- 
essary extension of the above loop is the addition of a 
quadrature loop to the existing configuration. Indeed, this 
is the case, as is illustrated in the quadriphase decision- 
feedback loop of Fig. 6-7. Letting <t>(t) represent the 
loop’s estimate of *J> (t) and 

r u (f ) = V~2 Ki cos $ (f ) 

r t (t) = -\T2K t sin$(f) (6-5) 

respectively, denote the upper and lower loop reference 
signals, then the upper and lower dynamic reference sig- 
nals €„ (t) and €; (t) are given by 

e«(t) — y-g- (V Pnd,(t) sin <f>(t) 

+ V Pd d 2 (t) cos 
+ V2N„[t^(t)]} 

( t ) = { - YPd di(t) cos <j>{t ) 

+ V Pod 2 (t) sin rf>(t) 

+ y[2 AUt,*(t)]} (6-6) 

where <f> (t) = <£ (t) — $ (t) and N u ft, <f> (t)] and N t [f, <f> (t)] 
are uncorrelated noise processes, which are modeled as 

N u [t, <j> (f)] = N c ( t ) cos <f> ( t ) — N s ( t ) sin <f> ( t ) 

Ni (t)] = N c ( t ) sin ( t ) + N s (t) cos <f> ( t ) (6-7) 

Proceeding with an analysis analogous to that presented 
for the suppressed carrier DAL, it can be shown that the 


signals z u ( t ) and z t ( t ) are given by 

2 “ w = K yy m Fi (p) ex p(~ p t ) 

X {ypPodt (t)di (t)sin<t>(t) 

+ ■\TPo d 2 (#) d, (t) cos 4> (t) 

+ V2 £ (t)N u [*,<£(*)]} 

Zi ^ = *Mp) ex p ( - p t ) 

X { — V Pod, ( t ) d 2 ( t ) cos <j> ( t ) 

+ V Pd d-, ( t ) d 2 ( t ) sin 4> (t) 

+ V2&(t)N l [t,<t,(t)]} ( 6 - 8 ) 

Following the same averaging argument for E (d,(t) di(t)} 
as previously given for the biphase DAL and the fact that 
E ( d , (t) d, (t)} = Ofor i=£j, (6-8) reduces to 

z„ ( t ) = F, (p) exp (- pT) 

X VPiTfl ~ 2 Fei [d> (*)] ] sin </>(#) 

+ \T2 d,(t)N u [#,^(t)]} 

Zl ® = F * (P) ex P ( ~ P T ) 

X {yfPo [1 - 2F E2 [<£(*)]] sin <£(#) 

+ 'f2&(t)N l [t,<t>m ( 6 - 9 ) 

It is clear that, in general, the transfer function factor 
exp (— pT), with p = /«, affects loop stability and reduces 
the signal acquisition or pull-in range. These phenomena 
have been discussed relative to the binary DAL, and thus, 
we shall not pursue them further here. Instead, we 
make the usual assumption that W L T < < 1, whereupon 
exp ( — ;<oT) is approximately unity for all <o within the 
equivalent bandwidth of the loop. Thus, from a steady- 
state performance standpoint, this factor has negligible 
effect. 
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Fig. 6-7. A quadriphase decision-feedback loop 


Letting z ( t ) — z u (t) + Z\ ( t ), then the stochastic integro- 
differential equation of the quadriphase decision-feedback 
loop becomes 

4>(t) = 6(t) 

(p) [1 - 2 P E1 [* (*)]] sin * ( t ) 

+ yf2Z(t)N u [t,4>(t))} 

~ yj F 2 (p) {VM1 - 2P e , [* (<)]] sin <t> (t) 

+ \[2 d 2 (t) Ni [t,4>(t)]} (6-10) 

A A 

where K u = K t K lm K v and K x - K,K 2m K v . 

The quadriphase decision-feedback loop of Fig. 6-7 can 
in principle be extended to N phases To see this, we first 
note that the polyphase signal in (6-2) can also be written 
in the binary representation 

N/2-1 

s(t) = ( sin Jf) V 2P » d « (0 sin ( 

k - o 

( 6 - 11 ) 


where the set {d k (t); k = 0,1, 2, , N/2 — 1) are again 

± 1 digital waveforms whose transitions occur at T -second 
intervals. Equation (6-11) represents a generalization of the 
qua'driphase signal given in (6-3) The N-phase decision- 
feedback loop then would have N / 2 arms corresponding 
to the N/2 decisions d k ( t ); k — 0, 1, 2, , N/2 — 1 on 

the binary signals d k (t), k = 0,1,2, • ,N / 2— 1. The 

output of the VCO would be phase-shifted by k-rr/N; 
k = 0, 1, 2, , N/2 — 1 to provide N/2 coherent refer- 

ence signals. As a result, the N/2 equivalent additive low- 
pass noise processes at the input of the summing junction 
are no longer mutually uncorrelated, and hence, for N > 4, 
the loop becomes less efficient. 

Another N-phase tracking loop which makes use of the 
principle of decision feedback and is equally efficient for 
all N is illustrated in Fig 6-8. Here, in each T-second 
interval, the decision 0 L on the transmitted phase symbol 
0 k — (2k + 1) 7 t/N is used to produce the decision-feedback 
signals. Under the same assumptions given for the pre- 
vious decision-feedback loop, the dynamic error signal 
becomes 

e(f) = K,K m { \[Ti> cos (0 k - 0 k ) sin <f, (t) 

+ yfPj, sin (6 k - fi k ) cos $ ( t ) 

+ cos $ k N„ [t,<f>(t)] + sin&N, [f,<£(t)]'} 

( 6 - 12 ) 



Fig. 6-8. An N-phase decision-feedback loop 
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where N u [t,4>(t )] and N, [*,<£(*)] are given by (6-7). 
Based upon the same averaging argument used previ- 
ously, the output of the loop filter in the tracking mode 
can be expressed in terms of the circular moments of 
- Ok, viz., 

z(t) = K,K m F(p) {VPdCOs( 0* - £&) sin <j> (t) 

+ sin (6 k - tk) cos 4> ( t ) 

+ cos 'O'kNu [t,4>(t)] + sin $ k Ni [t,4>(t)]} 

(6-13) 

This discrete random variable, d k ranges over the set 
of allowable values 2 /tt/N, j = 0, ±1, ±2, , ±N/2 — 1, 

N/2 with probabilities (Ref 6-1) 

a ( 2/V ) 

p > (<*>) = pr« b | A* - 8* = ~jfj 

= —J |^exp { — (u — YRdCOS</>) 2 } 

/ « tan [ (2j +Uir/lV] “I 

exp ( — (v — \TR d sin 4>) 2 } dv du 

utan [(2J-1 )tt/2V] J 

(6-14) 


where the prime on the summation denotes omission of 
the j = 0 term. 

We note that P 0 (<f>) is the conditional probability that 
the decision on B k is correct given the phase error <j>. Thus, 
letting 

N/2 

P E (*) = 1 - Po (4>) = ^' P, (*) (6-17) 

)=-N/ 2+1 

Eq. (6-16) may be expressed in the equivalent form 

N/ 2 

sin (0 k - $ k ) = "y ^ P, ( <t > ) sin ■— 

J--N/ 2+1 

cos (Sfc — &) = 1 — 2Pe (<#>) 

N/2 

+ P, W 1 + cos 

;=-iv /2 + i 

(6-18) 

Substituting (6-18) into (6-13) and recalling that 


where we have assumed that the loop phase error <f> ( t ) is 
essentially constant over several signaling intervals. Also, 
from the law of total probability. 


^2 p > w = 1 


(6-15) 


4>(t) = 0 (t) — Ky ~~~~ (6-19) 

the stochastic integro-differential equation of operation for 
the JV-phase decision-feedback loop of Fig. 6-8 becomes 


<f>(t) = 0 ( t ) - KF (p)i V Pd [1 - 2 P E (*)] sin * 


Thus, from (6-14) and (6-15), the circular moments of 
Ok ~ Ok can be expressed as 

V/2 

7T A 2/V 

sin (0 k - 0 k ) = y P, (<f>) sin -jj- 

J=-N/ 2+1 
N/2 

E ' 2/V 

Pj (</>) sin -jj- 

)=-N/ 2+1 


N/2 

+ "NTPo y ^ P, ( <j > ) sin <f> ( 1 + cos 

J = -N/2+1 
N/2 

+ V?D y ^ p, (4>) cos sin 

J=^N/2+l 

+ cos 'Ok N u (t, <f>) + sin 0 k N , (t, </>) 1 (6-20) 


) “v A „ , . 2/V 

cos (0 k - 0 k ) = y P, (<f>) cos -jj- 

J= -N/2+1 

y/g 

= Po fa) + P, fa) cos ^ 


; = -.v /2 + i 


(6-16) 


A 

where K = K a K m Kv. Recognizing from (6-14) that P, fa) = 
P -7 ( — <j>), the second and third terms of (6-20) are odd 
functions of 4> and as such contribute to the overall track- 
ing error characteristic. When N = 2, both of these terms 
are zero and (6-20) reduces to the stochastic integro- 
differential equation of operation for the biphase sup- 
pressed carrier DAL as given by (3-15). 
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In conclusion, we note that the theory presented in this 
section and the results given in Ref. 6-2 can be used to 
study the cycle slipping phenomenon of these two cir- 
cuits. This, however, shall not be done here Also, a dis- 
cussion of the data detection and noisy reference problem 
associated with the use of polyphase signaling is pre- 
sented in Ref. 6-1. 

E. Single-Sideband Receivers 

The data-aided receiver may be used to track and de- 
modulate single-sideband suppressed carrier (SSB-SC) 
signals The mathematical form of such a signal is given by 

s (t) = V Pod ( t ) cos w c t — ypP~ D d ( t ) sin u> c t (6-21) 

where d ( t ) is the Hilbert transform of d ( t ) (See Ref 6-3 
for a general discussion of the SSB-SC signal.) 

There are two ways by which the DAR can be realized 
The first involves choosing the input BPF (see Fig. 2-5) so 
that it is essentially centered on the sideband If an ideal 
rectangular filter is postulated, then the outputs of the 
loop phase detectors (Fig 3-2) can be written as 

e, =^T' (VPod(f) cos <j> — V~Po d (t) sin + n„(t)} 

( 6 - 22 ) 

e 2 = {V"Pod(t) sin <f, - y[P~ D d (t) cos <£ - n 2s (f)} 

(6-23) 

where n u (t) and n 2s (<) are the low-pass equivalents of the 
single-sideband noise formed by the input BPF operating 
on n, (t) (see Eq. 3-2), and are independent and gaussian 
with single-sided spectral densities N 0 /2 W/Hz. 

When <f> = 0, Eq (6-22) is essentially the same for the 
basic DAL discussed in Section III We see, however, that 
e 2 now contains an extra term, — V P D d(t), The presence 
of this term has no direct effect on the loop error function 
since E {d(t)} = 0 In fact, d(t) is a series of identical— 
except for sign— odd functions, each of which is centered 
on the respective symbols comprising d(t) Even when 
<j>^0, the presence of d in either c, or e 2 will have no 
effect since d, when averaged over one symbol period 
(excepting large symbol timing errors), is zero. Thus, the 
DAL will operate as if it were tracking a double-sideband 
suppressed earner signal. 


What counts in terms of the receiver performance is the 
exact variance of n ls and n 2g Ideally, as already stated, 
these noise components will have noise spectral density 
N 0 / 2. However, with practical bandpass filters, it is impos- 
sible to filter out all of the noise on the side of the carrier 
frequency for which the sideband has been suppressed 
Thus, the vanance of n ls and n 2s will be somewhat larger 
than ideally possible. For two-pole bandpass filters, the 
SNR degradation due to this effect has been measured to 
be about 0.7 dB The P*-related performance is shown in 
Fig. 6-9, together with the double-sideband performance 
taken from Fig. 4-19 

The second method by which the SSB-SC DAL receiver 
may be implemented is shown in Fig 6-10 The blocks 
labeled HT perform the Hilbert transform operation. For 
this configuration, there are no special requirements on 
the input BPF. Now, when <j> ( t ) is slowly varying with 
respect to d ( t ), it can be shown that 

e, + e" 2 = V 2 P D d (t) cos <j> -f m ( t ) — n 2 (#) (6-24) 

€2 + = V 2 Pod ( t ) sin (f> + n , ( t ) — n 2 ( t ) (6-25) 



Fig. 6-9. Single-sideband suppressed carrier DAR performance 
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Fig. 6-10. An alternative SSB-SC data-aided receiver 


where fh — n 2 and n t — n 2 are uneorrelated, gaussian, and 
have single-sided noise spectral densities 2N 0 W/Hz. This 
receiver is ideally optimum, and in practice will perform 
according to how well the Hilbert transform operations 
can be implemented. 

F. Subsymbol Integration and Decision for 
Very-Low-Data-Rate Applications 

In Section III-D, it was shown that when the symbol 
rate RR S approaches the loop noise bandwidth W L , the 
DAL will become unstable. From Eq. (3-60), we require 
that RR S >4.3W t if the loop is to be stable. Thus, if W,, 
is specified by, say, the need to track signal doppler rates, 
then there is a lower bound set on the symbol rate. For 
example, if W L = 18 Hz and RR„ = fR b = 16 bps (an un- 
coded bit rate), clearly the loop would be unstable. 

For the situation in which the symbol (or bit) error 
probability is nominally small, there is a solution to the 
problem. For the bandwidth and bit-rate numbers already 
cited, suppose that, in addition, P\ = 10~ 5 . These numbers 
are typical of a command situation. Consider now that 
each bit is subdivided into n equal-lengthy sub-bits, and 
that the DAL is implemented so that its a decisions are 
made on the basis of the sub-bits rather than the whole 
bit (Fig 6-11). Then, as far as the loop is concerned, it is 
operating at an effective symbol rate of nRR h , so that the 
stability criterion of (3-60) becomes 

nRR b > 4.3 W L (6-26) 

The cost is that the error probability which determines 
the efficiency loss L s (see Sections II and III) is much 
larger. However, if P\ is very small to begin with, this is 
a tolerable situation. Let us continue the example with 


the above command system numbers. Solving for n using 
Eq. (6-26) gives n > 4.8. Let us select n = 8, since this then 
gives a considerable stability margin. The effective symbol 
error probability for the loop then becomes 6 X 10~ 2 If 
we use the L a = 1 — 2 P* relationship, then L R = 0.6 dB 
and the loop SNR, p, is 8.8 dB, an acceptable operating 
point. The bit decision itself is, of course, still based on a 
whole bit, so that the bit error probability remains 10~ 5 . 

G. Two-Way Systems 

The application of DALs to two-way communication 
systems necessitates consideration of three possible 
modulation-demodulation alternatives 

(1) Suppressed carrier uplink and discrete carrier down- 
link— the spacecraft receiver contains a DAL and the 
ground receiver contains a discrete carrier-tracking 
loop. 

(2) Discrete carrier uplink and suppressed carrier 
downlink— the spacecraft receiver contains a dis- 
crete carrier-tracking loop and the ground receiver 
contains a DAL 

(3) Suppressed carrier uplink and suppressed carrier 
downlink— the spacecraft and ground receivers both 
contain DALs. 

Since Fig. 6-1 shows the necessary modifications of an 
existing spacecraft transponder receiver to accommodate 
a DAL, let us begin by considering the first of the above 
three alternatives Also, we shall assume two-channel 
configurations for both the up- and downlinks, wherein 
the transmitted signal is modeled as in Eq. (6-1). 

A block diagram of a two-way communication system 
based upon the above considerations is given in Fig. 6-12. 
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The ground transmitter emits the phase-modulated wave- 
form 

( t ) = V 2Pd„ (ill ( t ) COS a c t + 

+ Y 2 Po n d 2 1 (#)5 21 (t) sin < 0 C t (6-27) 



Fig. 6-11. Sub-bit decision OAR 



GROUND SYSTEM 


We observe m the spacecraft the phase-shifted, noise- 
corrupted waveform 

!/i (t) = V 2Pxi n dn (t) cos (oct + 0,) 

+ v 2 Pd 21 d 2 i (f)5 2 i (t) sin (<a c f + 0,) + n, (t) 

( 6 - 28 ) 

where n, (t) is white gaussian noise of single-sided spectral 
density N, u W/Hz Since a carrier component is not 
present in the received signal, the spacecraft employs a 
DAL to reconstruct a pair of reference signals r 28 (f) and 
^ 2 c (t) for demodulating dn (t) and d 2 i (t)5 21 (£) oflF the ear- 
ner The carrier reference signal r 28 (t) is also supplied to 
a subcarrier regenerating DAL (see Fig 6-5) for demodu- 
lating do, ( t ) off its subcarrier. 

If we model r,,(t) and r. c (f) as 

r 28 ( t ) = \T2sin [<u c * + $i (t)] 

( 6 - 29 ) 

T ic ( t ) = V~2" COS + 0, (f)] 

where 0i (£) is the DAL’s estimate of 9,, then on the down- 
link, we transmit* 

s 2 ( t ) = V 2Pr 2 cos [Ca> c t + 9„ u d 12 (t)Si 2 ( t ) 

+ O m „cU 2 (t) 5 22 ( t ) + G9 ; (t)] (6-30) 

‘Whether the standard or interplexed two-channel modulation for- 
mat is used here is of no consequence 



SR AC EC RAF! SYSTEM 


Fig. 6-12. Two-way communication system block diagram 
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where G is the transponder gain, whose value depends on 
the frequency band assignments of the uplink and down- 
link, e.g., S/S, S/X, X/X. At the ground station, we receive 
the phase-shifted, noise-corrupted waveform 


with W L 2 denoting the two-sided bandwidth of the 
ground receiver’s PLL. The necessary modifications of 
(6-34) for a second-order loop follow those given in Sec- 
tion III-C5 for the DAL. 


y 2 (t) — V 2P T 2 cos [ Goiot + @m 12 d, 2 (t) 5 12 (0 

+ Om 22 d 22 (t) 5 22 ( t ) + GO i (t) + 0 2 ] + n 2 (t) 

(6-31) 

where n, (t) is again white gaussian noise of single-sided 
spectral density N 02 W/Hz and is independent of n, (t) 
Separating the signal component of y 2 (<) into its spectral 
components, we get 

y , ( t ) = V2Pr 2 cos [G*, c t + GO, (t) + 0 2 ] 

- V 2Pd 12 di 2 (t) 5 12 (*) Sin [Ga> c t + GO, ( t ) + 0 2 ] 
— V 2 Po 22 d 22 (f)S 22 (t) sin [G<o c t + G0 t (f) + 0 2 ] 
— V 2P i 2 d 12 (t) d 22 (t) 5 12 (t) 5 22 (t) 

X cos [G Wc f + G& (f) + 0 2 ] + n 2 (t) (6-32) 

The ground station receiver tracks y 2 (t) with a phase- 
locked loop (PLL) and produces the reference signal 

r, (f) = V^sin [G«o c t + 0 2 (t)] (6-33) 

where 0 2 ( t ) is the PLL’s estimate of GO, (t) + 0 2 . 

The phase error probability density functions (p d.f.’s) 
in the spacecraft and ground receivers’ tracking loops are 
determined as follows Let 


<f>, ( t ) — 0, (t) 0, ( t ) 

and 

$2 (t) = 0 2 (t) 0 2 (f) 

be the loop phase errors in the spacecraft and ground 
receiver tracking loops, respectively, and let <j > , and </>_> 
be the modulo-2ir reduced versions of these loop phase 
errors. Then, p (<j>,) is given by (3-16) together with (3-17), 
where P D and R s are now replaced by P D „ and R, u = 
P ri ,T,,/N ul , with T„ denoting the reciprocal of the data 
rate of d„ (t). The p.d.f. p is given (assuming a first- 
order loop for simplicity) by the Tikhonov distribution 


where 


, , , exp (p 2 cos <f> 2 ) 
V \*P-I 0_I / 


2P C2 

pI ~ n 02 w lj 


(6-34) 


One problem which is becoming of increasing impor- 
tance in the study of two-way systems (particularly m the 
case of S/X-band transmission, where the transponder 
gain G is on the order of four) is the uplink command 
modulation feedthrough, which causes suppression of 
the downlink carrier. This problem has been investigated 
in Refs 6-4 and 6-5 with respect to its effect on Mariner 
transponders operating with two-channel command mod- 
ulation such as on the Mariner Venus Mercury 1973 
mission and on the Mariner Jupiter Saturn 1977 mission, 
which will use a single-channel type of command modu- 
lation In both of these cases, the spacecraft receiver 
employs a PLL for tracking the carrier, thus, as shown 
m Ref. 6-4, knowledge of the closed-loop transfer 
function of the PLL is sufficient to compute the down- 
link carrier suppression caused by the modulation com- 
ponents which are tracked by the PLL 

In the two-way system proposed here, the spacecraft 
receiver employs a DAL, for which the calculation of the 
downlink carrier suppression due to command modula- 
tion feedthrough requires a different technique. As dis- 
cussed in Section III-B, the input to the VCO is z ( t ), as 
given by (3-4), and which can be rewritten as 

z ( t ) = K 2 K 2m F (p) exp ( -pT„) V Pd„ { [d ( t ) d ( t ) 

- erf (V R«„ cos <)>,)] sin </>,} + K 2 K 2m F (p) 

X exp ( — pT«) {VTP^ erf (Vr^ cos <f> ,) sin </>, 

+ d(t)N 2 [t,4>,(t)]) (6-35) 

The second and third terms of (6-35) represent the signal 
and noise terms previously discussed in deriving the sto- 
chastic integro-differential equation of the DAL. The first 
term in (6-35) represents a self-noise due to the digital 
nature of the data detection process and the possibility 
of errors m the data caused by the additive noise and loop 
phase error Heretofore, this self-noise term has been 
neglected It behooves us here to investigate the effect of 
its presence in the phase of the reference signal produced 
by the VCO, since any modulation present here appears 
(multiplied by the factor G) on the downlink carrier 

Based upon arguments given in Section III-B, the self- 
noise process 


n s ( t ) = d(t)d ( t ) — erf (V R Sll cos <j>,) (6-36) 
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has zero mean and variance (conditioned on fa) given by 


ol, (fa) = E n { n 2 ( t )} = 1 - erf 2 ( cos fa) (6-37) 

(Note that as R e -> oo, -> 0.) Furthermore, since n t (t) 
is in essence a random pulse train, which in each sym- 
bol interval takes on value 1 — erf (V cos fa with 
probability 


1 + erf (V R Hl cos fa) 
2 


and value — 1 — erf(V cos $,) with probability 

l-erf(VX7 1 cos fa) 

2 


its power spectral density can be calculated based upon 
the techniques presented in Ref. 6-6. Denoting this density 
(conditioned on <£,) by S„ K (/; fa), we obtain 


(f, <frl) — rp 2 

n ~ -oo 

-f-(l -p) s 2 ( T7 , 01 )| 

+ rp p (1 p) | Si (/; fa) S 2 ( f ; <!>!) | (6-38) 

where S (x) is the Dirac delta function and 

P == [1 + erf (V R«„ cos fa)] 

Si (/, </>i) = [1 - erf (V R 8ll cos ^> t )] T s 
X exp ( - 

S 2 (f, fa) *= [-1 - erf (V RauCos<^i)] T e 

X«p(-MT.)[=^] (6-39) 


Substituting (6-39) into (6-38) gives, upon simplification, 

S». (/;*>) = T. [1 - erf 2 (\^cos<M][^^J 

(6-40) 


Note that 



(/, fa)df = 1 ~ erf 2 (V R.ii cos ^) 


(6-41) 


which agrees with (6-37) The self-noise process is passed 
through the loop filter, with the transfer function given by 
(3-46) Since the bandwidth of this filter is quite narrow 
relative to that of the self-noise power spectral density 
S„ s (/; fa), we can model F (p) as a rectangular low-pass 
filter (LPF) of unit height and width 1/t, Then, at the 
input to the VCO, the self-noise appears as a smooth, 
continuous type of random modulation given by 

«» ( t ) = K 2 K 2m V”P Dji sin fa [F(p)n,(t — T,)] (6-42) 


with approximate mean-squared value (conditioned on <f > ,) 

•*'.(*1 ) = £{«;*(*)} 

= (K 2 Kom V Fn,, sin </>, ) 2 £ 


1 - erf 2 (V R»„ cos </>i)~ 


(6-43) 


Recalling that the VCO is characterized by (3-14), then 
(#) appears as a frequency modulation on the downlink 
earner Since the rms value of n' K ( t ) is quite small, the fre- 
quency deviation of the downlink carrier is very small, 
and hence, the frequency spectrum of the modulated 
downlink carrier, due to the presence of the self-noise, is 
approximately given by 


if, fa) 


where 


7rW(0) 

[ttW (0)/2] 2 + (2irf — G<i>c) 2 


(6-44) 


w (0) = Ki Kk Kfi P Dii G 2 sin 2 fa [1 - erf 2 (VR^ cos fa)] 

(6-45) 

The spectrum in (6-44), which results whether the 
modulating self-noise n' s (t) is gaussian or not, has the bell- 
like shape of a Witch of Agnesi centered on the downlink 
carrier frequency Gw c , with a width ttW (0) between 3-dR 
points. Note from (6-44) that the spectrum of the down- 
link modulated carrier does not possess a discrete carrier 
component, but instead contains a continuum of values 
in a narrow band around the downlink carrier frequency. 
This phenomenon, namely, the absence of a discrete car- 
rier component in the modulated carrier spectrum, is not 
unlike the effect on the downlink carrier produced by the 
spacecraft transmitter instability itself, i e. the downlink 
carrier is never a pure line spectrum. However, as long as 
the bandwidth of the modulated downlink carrier around 
Gto c is small (relative to the bandwidth of the ground 
receiver’s tracking loop W L2 ), then the effective carrier 
suppression produced by the self-noise of the spacecraft 
tracking loop will be small In fact, we may define the 
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effective carrier suppression on the downlink due to the 
self-noise n’ s (t) by 


vo(4> i) 


J 

f\(f,4>)df-zl 

-00 J 1 

s. (/,*)# 

GU> c /2V + W L2 /2 


f °°S, (f; <l>) df 

I -co 


(6-46) 


Substituting (6-44) into (6-46) gives 



(6-47) 


The actual value of carrier suppression would be obtained 
by averaging (6-47) over the p.d.f of <f > i 

For large loop SNR, we can approximate sin <f> 1 by rf>, 
and cos <t> l by 1. Then, making use of (3-23) and (3-33), 
and noting that for small x, tan -1 x s x, we get 


A K 2 C 2 P 

Vc = VC (fa) = 1 - — ^ t 1 ~ erf2 (V R *„)] 

^1_ KGF £u _ [ i_erf 2(V ^-)] 

TTpl Wto 


K 2 G 2 N„ 


r 1 — erf 2 (V R.J 1 

2-n 


L erf 2 (\R,J J 


(6-48) 


where K = K,K. m K v , and the overbar denotes statistical 
average with respect to the p.d.f. of fa. For small P‘, 
Eq. (6-48) can be further approximated by 

2/W, \ 

vc^l- ) K 2 G 2 N 0 P‘ e (6-49) 


As an example of how small the carrier suppression is, 
consider the Mariner Jupiter Saturn example cited in 
Ref. 6-4, with the minimum operating conditions of K ^ 
4 X 10 7 and P* = 10~ 5 , together with N 0 ~ —167 dBm/Hz, 
G = 880/221, and W L2 /W L j ^ 1. Substituting these 
numerical values into (6-49) gives 

1(! sl- 10- 8 (6-50) 

In contrast to the worst-case condition for the discrete 
carrier system (carrier suppression = 0 83), which occurs 
for strong signal conditions, the worst-case condition for 
the DAL receiver occurs at system threshold, and, as can 
be seen, is totally negligible. For strong signal conditions 
with the DAL, P\ -» 0 and therefore, j j c — > 1. This supports 
our earlier claim that the DAR receiver does not exhibit 
any significant feedthrough problems 

In conclusion, we remind the reader that we have con- 
sidered only one of the three modulation-demodulation 
alternatives cited in the beginning of this subsection. 
Extension of the foregoing results to the other two alterna- 
tives is straightforward. 
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VII. Summary and Conclusions 


In summary, we have presented an efficient means for 
coherently demodulating a suppressed carrier signal 
which has been modulated by digital data. The data may 
either be in the form of binary symbols (or bits), in which 
case they are transmitted as biphase modulation of the 
carrier, or they may be N-level information, in which case 
the carrier is polyphase-modulated. In either situation, the 
carrier -tracking loop, which reconstructs the demodula- 
tion reference signals from the received signal m additive 
white gaussian noise, employs the principle of decision- 
directed feedback; i e., the operation of the loop depends 
upon the data modulation detection decisions (see 
Table II-2). This dependence greatly enhances the carrier- 
tracking ability of a loop which is required to track a 
data-modulated suppressed carrier. Clearly, an increase 
in the efficiency of carrier tracking manifests itself as an 
improvement in the area of radio link-derived navigation, 
and in the performance of the data detection, since the 
reference signals supplied by the loop contain less noise 

From the theoretical and experimental results pre- 
sented m this report, we are able to conclude that, for a 
wide variety of applications, the concept of suppressed 
carrier transmission together with a data-aided tracking 
receiver is generally more efficient than its discrete carrier 
counterpart which employs a PLL Despite this, we must 
proceed with some caution since a suppressed carrier 
data-aided tracking receiver does possess certain anoma- 
lies relative to the tracking receivers which are employed 
in present communication system designs. Furthermore, 
the operational experience with discrete carrier systems 
cannot be discounted. We shall review these anomalies 
here, comment on their significance, and discuss how their 
potential detrimental effect may be reduced by the tech- 
niques presented in the text. 

In the first case, suppressed carrier binary transmis- 
sion (as opposed to discrete carrier binary transmission) is 
inherently accompanied by a 180-deg phase ambiguity in 
the true system lock point This fact is true irrespective of 
whether a DAL or any other type of suppressed carrier 
tracking loop is employed. The significance of this am- 
biguity is with respect to its effect on the data estimates 
emanating from the data detector output With proba- 
bility 1/2, the data estimates may be sign inverted relative 
to their transmitted sense This problem can be resolved 
by transmitting a known synchronizing sequence (and 


recognizing it), or by differentially encoding the data 
(which results in a small penalty of about a factor of 2 in 
average symbol error probability performance). 

Secondly, the data-aided tracking loop incorporates a 
delay element which can potentially cause instability 
problems and which reduces the acquisition range of the 
loop. This delay also increases the loop bandwidth and 
natural frequency of operation relative to that which 
would exist if the delay were zero. These problems are 
discussed in Sections III-D and III-E From the results 
given there, we conclude that the problems introduced 
by the presence of the delay element are significant only 
for symbol rates comparable in magnitude to the loop 
bandwidth Furthermore, for such low data-rate applica- 
tions, and sufficiently small symbol error probability, one 
can employ subsymbol integration (Section VII-F), which, 
insofar as the tracking loop is concerned, increases the 
effective symbol rate by a factor of n, where n is the num- 
ber of subsymbols per symbol. This minimizes the delay 
and therefore its related problems 

Another seemingly potential problem is the effect of 
the modulation feedthrough in two-way links employing 
a DAL, caused by the loop’s dependence on the data 
modulation As it turns out, the magnitude of this feed- 
through is essentially negligible relative to the values 
characteristic of a discrete carrier system (Section VII-D). 
Thus, this heretofore unknown effect actually is an advan- 
tage in favor of the DAL. 

We now summarize some of the salient advantages of 
data-aided receivers (DARs) To begin with, in the DAR, 
no subcarrier is required for single-channel modulation 
since the data may be modulated directly onto the carrier. 
This results in a considerable savings in receiver hardware 
because a subcarrier tracking loop is no longer necessary 
Furthermore, since the DAL requires symbol (or bit) syn- 
chronization, integration, and decision, this processing, 
which has heretofore been separate or apart from the 
carrier tracking function, now becomes an integral part of 
the receiver 

The operating point of the DAL (as measured by loop 
SNR, etc.) does not depend on the structure of the data 
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symbols, and works equally well if long strings of ones or 
zeros, or a random mixture is transmitted. Nor does the 
DAL require any special acquisition prefix or sequence or 
comparative threshold decisions to achieve phase lock. 
(We have already noted the special problem of resolving 
the carrier phase ambiguity.) 

Since the proposed mechanization of the DAR involves 
sampled-data digital processing, data rate switching is 
easily accommodated. Sections IV and V deal extensively 


with implementation, including the important topics of 
automatic gain control and lock detection 

Finally, the section on complete system configurations 
(VI), gives considerable insight into how the basic DAL 
can be integrated into present receiver designs. In addi- 
tion, a number of common applications are discussed, and 
some suggestions are made concerning how the DAL may 
be used for bandwidth compressive modulations such as 
polyphase and single-sideband 
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Nomenclature 


A 

first level of quantizer output 

n a (t) 

self-noise process due to the digital na- 

d(t) 

data stream with ±1 symbols 


ture of the data detection process 

d(t) 

data detector’s estimate of d (t) 

No 

noise spectral density, W/Hz 

d(t) 

Hilbert transform of d (t) 

N, [*,*«], 

N 2 [*,*(*)] 

low-pass noise processes at the in-phase 
and quadrature phase detector outputs 

F(s) 

transfer function of loop filter 


of the DAL 

F, 

ratio of r 2 to tj 

P 

Heaviside operator 

G 

transponder turnaround ratio 

P(<t>) 

phase error probability density function 

G(<t>) 

difference between the loop S-curve and 

Pc 

probability of correctly determining 

its statistical mean with respect to p (<j>) 


loop’s lock condition 

G d 

ratio of r^ c to t]oc (data signal improve- 

Pc 

discrete carrier power, W 


ment of the suppressed carrier system to 
the discrete carrier system) 

Po 

data signal power, W 


Pt 

total transmitter power, W 

GtR 

ratio of p sc to p DC (tracking loop improve- 
ment of suppressed carrier-tracking loops 

PI 

system average symbol error probability 


to discrete carrier-tracking loops) 

r 

loop damping factor (equal to four times 

H(s) 

K 

closed-loop transfer function 

r(t) 

the square of the true loop damping) 
carrier reference signal produced by car- 

number of bits of quantization in ADC; 
product of K 2 , K 2m , and K v (total open- 

rier regeneration loop 


loop gain) 

R, 

data signal-to-noise ratio 

K (V) 

gain of controllable IF amplifier in AGC 

fRb 

system bit rate, bits/s 


loop 


data symbol rate, symbols/s 


a shape factor determined by transfer 

s(t) 

transmitted signal 


function of the bandpass filter in squar- 



ing and Costas loops 

S 

average rate of cycle slipping, slips/s 

Ky 

VCO gam 

Sn, {f, <f>l) 

power spectral density of n s (t) condi- 



tioned on <f>i 


In-phase and quadrature phase detector 
gains 

ef {<!>) 

loop S-curve (effective loop nonlinearity) 

K?, K 2 , 

average power of the in-phase and quad- 
rature reference signals provided by the 

5q(<£) 

a unit power square wave in <f> with 
period 2 tt 


DAL 

quantizer loss function for ADC 

5 (t) 

unit power square-wave subcarrier of 
angular frequency o>sc 

FaDC 

T d 

efficiency loss associated with the carrier 

arbitrary loop delay 

Lr 



regeneration process in suppressed car- 

T, 

reciprocal of £R a (data symbol period), s 


rier systems 

UoM 

loop potential function 

m a 

number of symbols of accumulation in 

V 

first threshold of quantizer input, also 


forming lock detector signal 


AGC voltage 

n, (t) 

additive input noise (narrowband gaus- 

V' 

threshold for testing normalized lock de- 


sian process) 


tector signal 
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Nomenclature (contd) 


Vi 

lock detector signal 

d m 

data signal modulation index 

V 9 

quadrature counterpart of AGC voltage 

00 

random phase component of 9 ( t ) 


used in lock detection 

m 

ratio of <*, at the design point to its value 

V? 

normalized lock detector signal 


at the operating point 

w, 

input IF filter bandwidth, Hz 

MV 

mean of AGC voltage 

w L 

two-sided loop bandwidth, Hz 

/4 

normalized mean of AGC voltage 

x(t) 

total received signal 

Uytfj (Tylt 

j /? 

mean and variance of V* 

y 

ratio of the design point loop bandwidth 

a 0. 

input signal-to-noise ratio to the ADC 

to the input IF filter bandwidth 

Padc 

z(t) 

input signal to VCO 

pDC 

loop signal-to-noise ratio for the discrete 
carrier-tracking loop 

a, 13 

parameters characterizing the phase 
error probability density function p ( <j > ) 

pi 

signal-to-noise ratio in the noise band- 
width of the input bandpass filter 

«i 

bandpass limiter suppression factor 


y (“/) 

synchronization boundary of loop 

pxc 

loop signal-to-noise ratio for the sup- 
pressed carrier-tracking loop 


limiter performance factor affecting loop 
signal-to-noise ratio 

<T% 

variance of the phase error 

a 

reciprocal of the single-sided loop band- 

a n* (0 1) 

variance of n s ( t ) conditioned on 0, 


width-symbol time product 

Tl, T 2 

loop filter time constants 

e 

epoch of transmitted data stream, s 

Yo 

design point of loop preceded by a band- 

€ 

symbol synchronizer’s estimate of € 


pass limiter 

e(t) 

phase detector output 

0 (0 

loop phase error 

Vr 

critical value of normalized loop delay 

*(*) 

total phase of received carrier 


determining absolute stability region 

$(#) 

carrier regeneration loop’s estimate of 

VC 

average of vc(<t>i) over the probability 


*(*) 


density function of <j> i 


steady-state phase error (static phase 

VC (01 ) 

effective carrier suppression on the down- 
link due to the self-noise process (condi- 
tioned on <f> i) 

0)C 

error) 

carrier angular frequency, rad/s 

Vdc 

Vsc 

efficiency of discrete carrier system 
efficiency of suppressed carrier system 

Of 

fundamental beat note in a Fourier ex- 
pansion of 0 ( t ) 

subcarrier angular frequency, rad/s 

<OSC 

e(t) 

phase shift added to transmitted signal 

by the channel characteristic 

Oo 

doppler frequency offset, rad/s 
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INDEX 


A 

Accumulator, 29, 32, 38 
Acquisition range, 24 
Acquisition 

DAL, 7, 11, 23-26 
Symbol synchronizer, 30-31 
Adder, binary, 38 
Amplifiers, RF, IF, 43, 48 
Analog mechanizations, 28 
Analog-to-digital converter ( ADC ) , 28, 29, 
35-39 

Anti-modulation device, 5 

Automatic gain control ( AGC ) , 8, 43-44 

B 

Bandpass limiter, 26-27 
Bandpass filter, 6, 50, 55 
Bandwidth 
Loop, 3, 10, 20 
Low-pass filter, 32, 34 
Radio frequency, 4, 5 
Sampling, 32-35 
Bandwidth savings, 5 
Bang-bang servo, 30 

c 

Carrier regenerating receiver, 2, 3 
Command, 3, 50 
Command detector, 7, 38 
Costas loop, 5-7 
Cycle slipping, 11, 16-17 

D 

Damping parameter, 18-20 

Data-aided loop (DAL), 4, 5, 7, 8, 9-27, 28 

Data-aided receiver (DAR), 2, 9, 28—41 

Data detector, 6, 10 

Data estimate, 6, 10, 28 

Data rate switching, 7, 28, 50-51 

Data signal improvement, 3 

Data stream, 3, 29-30 

Delay 

Effect on loop bandwidth, 20-23 
Effect on loop natural frequency, 20-23 
Effect on loop stability, 18-20 
Delay function, 8, 18, 31 
Delay function, mechanization, 31-32 
Delay lines, 31-32 
Design point, 26 
Demodulation, coherent, 6 
Detuning, 10 
Digital logic, 38 

Digital mechanizations, 28, 31-32 
Digital-to-analog converter (DAC), 28, 38, 
44 

Discrete carrier, 1, 3 
Doppler offset, 10, 24 
Downlink, 5, 56-60 


E 

Efficiency 

Analog-to-digital conversion, 37-38 
Complete DAL receiver, 41 
Discrete carrier signals, 3 
Loss, 3, 4, 6, 7 
Matched filter, 29 
Sampling, 32—35 

Suppressed earner receiver, 6-7, 10 
Suppressed earner signals, 3 
Symbol synchronizer, 31 
Two-channel signals, 4 
Envelope AGC, 8, 43-44 

F 

False lock, 23-26 
Filter, Butterworth, 6, 35 
Filter 

In-phase, 29 
Mid-phase, 29 
Filter, loop (see loop filter) 

Frequency offset 

Camer (see doppler offset) 

Timing, 31 

G 

Gain control, automatic (see automatic gain 
control) 

Gain, loop, 11, 53 
Gain, multiplier, 10 
Gain, VCO, 11 
Gaussian noise, 10 

H 

Hilbert transform, 55-56 
Hybrid DAL, 9, 51 

I 

Interplex, 4 

L 

Limit cycle, 23 

Limiter, bandpass (see bandpass limiter) 
Linear stability analysis (see loop stability, 
linear) 

Lock detector, 2, 8, 25, 44—46, 51 
Lock, false, 23-26 
Loop filter, 10 

First-order loop, 11 
Second-order loop, 17 
Loop nonlinearity, 12-13 
Loop stability, 11 
Linear, 18-20 
Nonlinear, 20 


M 

MPSK tracking, 50-55 
Matched filter, mechanization, 28-30 
Modifications of discrete carrier receiver, 
48-49 
Modulation 
Biphase, 3, 9 
Discrete camer, 3, 9 
Interplex, 4 
Phase, 3 
Polyphase, 51 
Quadnphase, 4, 51 
Single-channel, 3 
Single-sideband, 55 
Subcarrier, 3, 9 
Suppressed earner, 3 
Two-channel, 4, 50 
Modulation Feedthrough, 8, 58-60 
Multiplier gain, 10 

N 

Noise, bandpass, 10 
Noise, low-pass, 10, 28 
Noise, self-, 58-59 
Noise spectral density, 3, 10 
Nonlinear stability analysis (see loop stabil- 
ity, nonlinear) 

Nyquist plot, 18-20 
Nyquist rate, 29, 36 

P 

Performance tests, 41 
Phase ambiguity, 8, 13, 48 
Phase detector, 9, 25, 28 
Phase error 
Loop, 10 

Static (see static phase error) 
Steady-state, 12 
Phase-locked loop, 1, 9, 49, 56 
Phase noise, variance, 11, 13-15 
Potential function, 11-12, 17 
Power 

Data signal, 3 
Discrete carrier, 3 
Sideband, 3 
Transmitter, 3 
Power spectrum, 4, 50, 59 
Probability density function 
of phase error process, 11-13, 17 
of symbol synchronizer timing offset, 30 
Probability of detecting lock, 46 
Probability of symbol error, 10, 15, 41 

Q 

Quadriphase, 4 

Quadriphase tracking loops, 51-55 
Quantization loss, 35 
Quantizing characteristic, 36 
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INDEX (contd) 


R 

Ranging, 50-51 
Reference signals, 10 


s 

S-curve, 9, 12-13, 25 
Signal-to-noise ratio 
ADC loss, 37 
Loss in sampling, 34-38 
Square-law, 6 
Symbol, 10, 14 
Tracking loop, 3, 9 
Single-sideband DAL, 55-56 
Squaring loop, 4, 5-7 
Squaring loss, 6, 14 


Stability, loop (see loop stability) 

Static phase error, 12 

Stochastic integro-differential equation, 

9-11 

Subcarrier, 3, 4, 49-50 
Subcarrier demodulator assembly, 7 
Subcarrier tracking, 49-50 
Sub-symbol integration, 56 
Suppressed carrier, 1, 3 
Suppressed carrier receivers, 5-7 
Suppression factor, limiter, 26 
Symbol period, 4, 7, 10, 29 
Symbol rate, 4, 10 

Symbol synchronizer, mechanization, 28-31 
Synchronization boundary, 23 


T 

Telemetry, 2, 3 

Timing error, 28 

Tracking loop improvement, 4 

Transition detector, 29 

Transponder, 5, 56 

Turnaround ratio, 5, 57-58 

Two-way communication system, 56-60 

u 

Uplink, 5, 48, 56 

V 

Voltage-controlled oscillator (VCO), 6, 11, 
28 
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